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Chapter  1 


MICRO/MINIATURE  HEAT 
PIPES  AND  OPERATING 
LIMITATIONS 

1.1  Summary 

A  literature  review  related  to  miniature  and  micro  heat  pipes  is  given.  It  is  found 
that  longitudinal  groove  designs  are  crucial  to  increase  the  heat  transport  capacity  of 
miniature  heat  pipes.  In  addition  to  the  operational  limitations  of  conventional  heat 
pipes,  micro  heat  pipes  may  be  subject  to  the  vapor  continuum  limitation,  which 
may  prevent  micro  heat  pipes  from  operating  under  lower  working  temperatures.  An 
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analysis  of  the  capillary  limit  reveals  that  the  disjoining  pressure  may  play  a  role  in 
heat  transfer  for  the  micro  heat  pipe  and  may  increase  its  heat  transfer  capacity. 

1.2  Literature  Review 

Thermal  management  is  one  of  the  most  critical  technologies  in  electronic  product  de¬ 
velopment,  and  directly  influences  the  cost,  reliability,  and  performance  of  the  finished 
product.  As  the  number  of  circuits  on  a  computer  chip  increases,  it  becomes  more 
difficult  to  dissipate  the  heat  which  is  generated.  The  chip  heat  fluxes  in  mainframe 
computers,  for  example,  have  reached  60  W/cm2,  and  are  expected  to  exceed  100 
W/cm2  by  the  year  2000  (Jacobs  and  Hartnett,  1991).  In  addition  to  the  limitations 
on  the  maximum  chip  temperature,  there  may  be  further  requirements  on  the  level  of 
temperature  uniformity.  Because  of  the  high  heat  fluxes  and  temperature  uniformity 
considerations,  it  is  necessary  to  develop  new  methods  for  distributing  and  removing 
heat  from  modern  electronic  devices.  The  micro  heat  pipe  is  one  of  the  promising 
technologies  for  the  achievement  of  high  local  heat  removal  rates  and  uniform  tem¬ 
peratures  in  computer  chips.  Micro  heat  pipe  structures  can  be  fabricated  on  the 
surface  of  the  substrate  of  electronic  chips  using  the  same  technology  that  forms  the 
circuitry.  These  thermal  structures  can  be  an  integral  part  of  the  electronic  chip  and 
remove  heat  directly  from  the  area  where  maximum  dissipation  occurs  (Jacobs  and 
Hartnett,  1991).  Micro  heat  pipes  can  also  find  applications  in  cooling  the  leading 
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edge  of  stator  vanes  in  turbines  and  in  space  waste  heat  radiator  panels  by  increasing 
the  effective  thermal  conductivity  of  phase-change  materials. 

Cotter  (1984)  was  the  first  to  propose  the  micro  heat  pipe  concept  for  the 
cooling  of  electronic  devices.  The  micro  heat  pipe  was  defined  as  a  heat  pipe  in  which 
the  mean  curvature  of  the  vapor-liquid  interface  is  comparable  in  magnitude  to  the 
reciprocal  of  the  hydraulic  radius  of  the  total  flow  channel.  Typically,  micro  heat 
pipes  have  convex  but  cusped  cross  sections  (for  example,  a  polygon),  with  hydraulic 
diameters  in  the  range  of  10  to  500  /z m. 

Based  on  the  original  idea  proposed  by  Cotter  (1984),  a  silver  and  a  copper 
flat  micro  heat  pipe  with  dimensions  0.6  (W)  x  0.1  (T)  x  25  (L)  mm  were  manufactured 
by  the  ITOH  Research  &  Development  Corporation  (Japan).  A  heat  transport  rate 
of  0.3  W  was  reportedly  obtained,  which  corresponds  to  a  heat  flux  on  the  order  of 
1  W /cm2.  Flat  heat  pipes  with  dimensions  on  the  order  of  10  (W)  x  10  (T)  x  50  (L) 
mm  were  also  built  for  the  cooling  of  semiconductor  elements  (Fig.  1.1).  The  flat 
heat  pipe  can  be  installed  on  a  semiconductor  element  or  a  semiconductor  substrate 
with  the  flat  surface  of  the  heat  pipe  parallel  to  the  substrate.  On  the  inner  surface 
of  the  flat  wall  which  is  in  contact  with  the  semiconductor  surface,  small  grooves 
were  made  in  the  longitudinal  direction  to  provide  capillary  pumping.  The  heat  pipes 
manufactured  by  ITOH,  by  definition,  are  not  actually  micro  heat  pipes,  but  are 
in  the  miniature  or  semi-micro  range.  Since  most  microelectronic  components  are 
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Figure  1.1:  The  flat  heat  pipe  developed  by  ITOH  Research  and  Development  Cor¬ 
poration 
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shaped  as  thin  plates  with  flat  base  plates,  and  the  thermal  conductivity  of  silicon  is 
rather  high,  the  heat  transfer  resistance  in  the  transverse  direction  is  relatively  small. 
Therefore,  a  flat  heat  pipe  attached  on  the  component  surface  is  advantageous  in 
cooling  microelectronic  chips  with  moderate  or  high  heat  fluxes. 

Babin  et  al.  (1989)  tested  a  trapezoidal  heat  pipe  with  a  square  cross  section 
of  1  mm2  and  a  length  of  57  mm  manufactured  by  ITOH,  which  is  shown  in  Fig.  1.2. 
The  case  material  was  either  copper  or  silver  and  the  working  fluid  was  ultrapure 
water.  The  triangular  grooves  in  the  four  corners  of  the  heat  pipe  served  as  arteries 
for  the  return  of  liquid.  The  maximum  heat  transport  rate  was  about  0.5  W  at 
working  temperatures  of  60  to  70°C,  which  corresponds  to  heat  fluxes  of  about  1 
W/cm2  on  the  outer  surface  of  the  evaporator.  Wu  and  Peterson  (1990)  tested  a  flat 
heat  pipe  with  a  cross  section  of  1  mm  x  2  mm  and  a  length  of  60  mm  with  ultrapure 
deionized  water  as  the  working  fluid.  The  heat  pipe  was  specifically  designed  for  use 
in  the  thermal  control  of  ceramic  chip  carriers.  In  this  application,  the  heat  pipe  is 
fitted  securely  under  the  chip  and  is  attached  to  the  chip  carrier.  A  maximum  heat 
transport  rate  of  1.3  W  was  obtained  which  corresponds  to  a  heat  flux  of  5  W/cm2. 

Plesch  et  al.  (1991)  tested  two  types  of  miniature  fiat  heat  pipes  for  heat 
removal  from  microelectronic  circuits,  both  having  overall  dimensions  of  7  (W)  x  2 
(T)  x  120  (L)  mm.  The  first  type  of  heat  pipe  (Fig.  1.3(a))  has  two  arteries  which  are 
formed  by  the  spars  of  an  inserted  frame  and  the  wall  of  the  casing.  Transverse  grooves 
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were  provided  to  allow  transverse  transport  of  liquid  working  fluid  in  the  heating  and 
cooling  zones.  The  second  type  (Fig.  1.3(b)),  however,  has  longitudinal  grooves  over 
its  entire  length.  For  the  first  type,  a  maximum  heat  transport  rate  of  about  16.5  W 
was  obtained  in  a  horizontal  arrangement  with  a  cooling  water  temperature  of  67°C. 
The  temperature  drop  over  the  heat  pipe  length  in  this  case  was  about  25°C.  For 
the  second  type,  the  maximum  heat  transport  rate  obtained  was  nearly  70  W,  which 
corresponds  to  a  heat  flux  of  about  35  W/cm2  in  a  horizontal  arrangement  with  a 
cooling  water  temperature  of  42°C.  The  temperature  drop  over  the  heat  pipe  in  this 
case  was  about  35° C.  For  the  vertical  arrangement  of  the  heat  pipe  with  the  cooling 
zone  on  top,  the  heat  transport  capacity  was  even  higher,  where  the  maximum  heat 
flux  reached  60  W/cm2.  Compared  to  the  miniature  heat  pipes  built  by  ITOH,  which 
were  tested  by  Babin  et  al.  (1989),  Wu  and  Peterson  (1990),  Itoh  and  Polasek  (1990a; 
1990b),  the  maximum  heat  fluxes  obtained  by  Plesch  et  al.  (1991)  were  much  higher. 
This  difference  may  be  due  to  the  different  inner  surface  configurations  of  the  heat 
pipes,  which  provide  capillary  forces  for  liquid  to  flow  back  to  the  evaporator.  For  the 
heat  pipe  tested  by  Plesch  et  al.  (1991),  many  fine  axial  grooves  were  manufactured 
on  the  inner  wall  of  the  surface  which  contacts  the  electronic  device.  For  the  heat 
pipes  by  ITOH,  however,  only  four  corners  provide  capillary  pumping.  Since  the 
circumferential  length  of  the  heat  pipe  cross  section  is  relatively  large,  the  capillary 
force  provided  by  the  four  corners  is  obviously  not  large  enough  to  pump  a  sufficient 
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amount  of  liquid  back  to  the  evaporator.  Further  investigations  have  been  completed 
recently  in  the  area  of  miniature  heat  pipes.  These  include  Kojima  et  al.  (1992), 
Lee  et  al.  (1992),  Chen  et  al.  (1992),  Zhou  et  al.  (1992),  Li  et  al.  (1992)  and 
Zhang  et  al.  (1992).  Kojima  et  al.  (1992)  developed  and  investigated  a  new  cooling 
system  with  miniature  heat  pipes,  which  collects  and  removes  heat  from  several  large- 
scale  integrated  circuits  (LSIs).  The  diameter  of  the  heat  pipes  studied  was  3  mm. 
The  cooling  system  was  mounted  on  every  surface  of  the  LSIs  and  cooling  fins  were 
connected  to  the  heat  pipes.  The  performance  of  the  cooling  system  with  miniature 
heat  pipes  was  several  times  higher  than  the  normal  mounting  style. 

Lee  et  al.  (1992)  reported  a  visual  and  experimental  study  on  the  two- 
phase  flow  and  heat  transfer  characteristics  of  small  two-phase  closed  thermosyphons 
with  wire  inserts.  The  inside  diameters  of  the  tested  thermosyphons  were  1.0  ~  3.6 
mm,  and  the  wire  diameters  were  0.4  ~  1.5  mm.  They  concluded  that  for  a  two- 
phase  closed  thermosyphon  made  of  a  very  small  diameter  tube  (ID  <  2.3  mm),  it  is 
necessary  to  have  wire  inserts  to  achieve  normal  operation  as  a  heat  transfer  element. 
Also,  the  heat  transfer  characteristics  of  such  two-phase  closed  thermosyphons  with 
wire  inserts  were  very  similar  to  those  of  two-phase  closed  thermosyphons  made  of 
large  diameter  tubes  without  wire  inserts. 

Chen  et  al.  (1992)  conducted  visualization  experiments  for  a  glass  heat  pipe. 
The  outside  diameter  of  the  glass  capillary  is  2  mm,  the  wall  thickness  is  0.25  to  0.3 
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mm,  and  the  total  tube  length  was  113  to  118  mm.  Four  basic  two-phase  flow  patterns 
were  observed:  (a)  liquid  projectile  flow,  (b)  slug  flow,  (c)  suspended  liquid  column, 
and  (d)  annular  flow.  The  concept  of  the  two-phase  flow  in  porous  media  employing 
modified  relative  permeabilities  was  used  to  calculate  the  maximum  performance  of 
a  micro  heat  pipe.  The  calculated  maximum  performance  for  vertical  orientation  is 
about  three  times  larger  than  that  for  horizontal  orientation. 

Zhou  et  al.  (1992)  experimentally  studied  the  operational  characteristics 
of  a  copper/ acetone  miniature  heat  pipe.  They  found  that  the  maximum  operating 
power  of  the  heat  pipe  was  limited  by  the  capillary  limit,  which  was  independent  of 
the  cooling  air  velocity.  However,  the  cooling  air  temperature  had  an  important  effect 
on  the  maximum  power. 

Li  et  al.  (1992)  experimentally  studied  cylindrical  miniature  heat  pipes 
with  inner  diameters  of  1.2  mm  and  2.5  mm,  and  lengths  of  100  mm  and  120  mm. 
They  found  that  the  effects  of  the  amount  of  the  working  fluid  on  the  heat  transfer 
performance  for  small  heat  pipes  were  larger  than  those  for  conventional  heat  pipes, 
and  the  influence  of  entrainment  on  the  capillary  limit  is  more  important  for  a  small 
diameter  heat  pipe  than  a  larger  one. 

The  minimum  cross  sectional  dimension  of  the  above  miniature  heat  pipes 
was  on  the  order  of  1  mm,  and  the  heat  pipes  were  attached  to  the  electronic  com¬ 
ponents  to  be  cooled.  Since  the  minimum  dimension  of  the  heat  pipes  is  not  very 
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small,  the  performance  characteristics  should  not  deviate  significantly  from  those  of 
conventional  heat  pipes.  Recently,  attempts  were  being  made  to  etch  micro  heat  pipes 
directly  into  the  silicon  and  use  them  as  thermal  spreaders  (Peterson,  1990;  Gerner 
et  ah,  1992).  These  heat  pipes  are  expected  to  have  hydraulic  diameters  on  the  order 
of  10  fi m,  which  are  classified  as  actual  or  “true”  micro  heat  pipes.  As  the  size  of 
heat  pipes  diminishes,  some  performance  characteristics  are  expected  to  be  different 
from  those  of  conventional  heat  pipes. 

Recently  Khrustalev  and  Faghri,  (1993)  developed  a  detailed  mathematical 
model  in  which  the  heat  and  mass  transfer  processes  in  a  micro  heat  pipe  (MHP) 
are  examined.  This  model  describes  the  distribution  of  the  liquid  in  an  MHP  and  its 
thermal  characteristics  depending  upon  the  liquid  charge  and  the  applied  heat  load. 
The  liquid  flow  in  the  triangular-shaped  corners  of  an  MHP  with  polygonal  cross 
section  is  considered  by  accounting  for  the  variation  of  the  curvature  of  the  free  liquid 
surface  and  the  interfacial  shear  stresses  due  to  liquid- vapor  frictional  interaction.  The 
predicted  results  obtained  by  Khrustalev  and  Faghri  (1993)  are  compared  to  existing 
experimental  data.  The  importance  of  the  liquid  fill,  minimum  wetting  contact  angle, 
and  the  shear  stresses  at  the  liquid- vapor  interface  in  predicting  the  maximum  heat 
transfer  capacity  and  thermal  resistance  of  the  MHP  is  also  demonstrated. 
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1.3  Vapor  Continuum  Limitation 


As  the  size  of  the  micro  heat  pipe  decreases,  the  vapor  in  the  heat  pipe  core  will  lose 
its  continuum  characteristics.  The  continuum  criterion  is  usually  expressed  in  terms 
of  the  Knudsen  number  (Cao  and  Faghri,  1992a) 


<0.01  continuum  vapor  flow 
>0.01  rarefied  or  free  molecular  flow 


(l.l) 


where  A  is  the  mean  free  path  of  the  vapor,  and  D  is  the  minimum  dimension  of  the 
vapor  flow  passage.  For  a  circular  vapor  space,  D  is  the  vapor  core  diameter.  The 
heat  transport  capability  of  a  heat  pipe  operating  under  the  rarefied  or  free  molecular 
vapor  condition  is  very  limited,  and  a  large  temperature  gradient  exists  along  the  heat 
pipe  length.  As  a  result,  the  micro  heat  pipe  will  lose  its  advantages  as  an  effective 
cooling  method  for  microelectronic  devices.  Therefore,  this  heat  transport  limitation 
is  referred  to  as  the  vapor  continuum  limitation  in  this  paper. 

The  mean  free  path  given  by  Tien  and  Lienhard  (1979),  based  on  the  kinetic 
theory  of  dilute  gases,  is 


A  = 


1.051  nT 


V2: 


7T  azp 


(1.2) 


where  a  is  the  collision  diameter,  and  k  is  the  Boltzmann  constant. 

By  combining  eqns.  (1.1)  and  (1.2)  and  applying  the  equation  of  state 
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p  —  pRT 


(1,3) 


for  pressure  p,  one  can  obtain  the  transition  density  from  continuum  to  rarefied  or 
free  molecular  vapor  flow  as 


P  = 


1.051  k 


(1.4) 


\/27r  a2RD  Kn 

Assuming  that  the  vapor  is  in  the  saturation  state,  the  transition  vapor 
temperature  corresponding  to  the  transition  density  can  be  obtained  by  using  the 
Clausius-Clapeyron  equation  combined  with  the  equation  of  state 


rp  Prf 

Ttr  =  —5  exp 
pR 


hfg 
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(1,5) 


where  prj  and  Trf  are  the  reference  (saturation)  pressure  and  temperature,  and  the 
vapor  density  p  is  given  by  eqn.  (1.4). 

Equation  (1.5)  can  be  rewritten  as  (Cao  and  Faghri,  1992a) 


In (T„,  R/Ptj)  +  (y-  -  jb)  =  0  (1-6) 

and  solved  iteratively  for  Ttr  using  the  Newton-Raphson/secant  method  (Pletcher  et 
ah, 1988). 

Figure  1.4  shows  the  transition  temperature  for  a  sodium  heat  pipe  as  a 
function  of  minimum  vapor  space  dimension.  When  the  minimum  vapor  space  di- 
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mension  reduces  to  less  than  1  mm,  the  transition  temperature  rises  sharply.  A  higher 
transition  temperature  will  make  the  startup  of  the  heat  pipe  from  the  frozen  state 
more  difficult.  However,  even  at  a  minimum  vapor  space  of  50  /z m,  the  transition 
temperature  is  nearly  1100  K,  which  is  in  the  normal  operating  temperature  range  of 
800  to  1400  K.  This  means  that  once  the  heat  pipe  startup  is  completed,  the  micro 
heat  pipe  will  be  operational  in  the  upper  operating  temperature  range. 

Figure  1.5  shows  the  transition  temperatures  for  water  and  methanol  as 
heat  pipe  working  fluids.  As  the  minimum  vapor  space  dimension  approaches  the 
micro  heat  pipe  range  ( D  <  100  y^m),  the  transition  temperature  for  water  rises  to 
more  than  50°C.  In  this  case,  a  micro  heat  pipe  with  water  as  the  working  fluid  is 
apparently  not  suitable  for  cooling  electronic  devices,  which  work  at  temperatures 
preferably  lower  than  50° C.  Even  if  the  working  temperature  of  the  micro  heat  pipe 
is  higher  than  the  transition  temperature,  when  a  electronic  device  starts  up  from 
the  room  temperature,  the  heat  pipe  will  temporarily  experience  the  rarefied  or  free 
molecular  vapor  flow  stage,  in  which  the  temperature  gradient  in  the  axial  direction 
is  substantially  higher  than  that  of  the  normal  working  condition.  This  continuum 
limitation  is  less  severe  for  methanol  as  the  working  fluid.  At  a  minimum  vapor 
space  D  =  50  //m,  the  transition  temperature  is  only  about  20° C.  However,  the  vapor 
pressure  of  methanol  is  higher  than  that  of  water.  At  a  working  temperature  of 
about  80°C,  for  example,  the  vapor  pressure  of  methanol  is  about  1.72  x  105  Pa,  in 


15 


Figure  1.5:  Transition  temperatures  for  water  and  methanol  as  heat  pipe  working 
fluids  (Cao  and  Faghri,  1994) 
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comparison  with  the  vapor  pressure  of  about  0.47  x  105  Pa  for  water.  A  working 
vapor  pressure  higher  than  the  ambient  pressure  may  result  in  additional  structural 
requirements  for  the  silicon  chip. 

Boiling  Limitation 

The  critical  temperature  drop  across  the  heat  pipe  wick  for  the  boiling  limit 
is  (Cao  and  Faghri,  1992b) 


AT, 


cut 


Ts-Tv  = 


R  TVTS 
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In 


fa 


l  +  W  J_  _  1  \  {  Pv 
Pv  \  n>  Rm  J  PvHPe. 


(1.7) 


where  Ts  is  the  wall-liquid  interface  temperature,  Tv  is  the  heat  pipe  vapor 
temperature,  r&  is  the  effective  radius  of  small  vapor  or  noncondensible  gas  bubbles 
trapped  at  the  inner  wall,  and  Rm  is  the  radius  of  the  liquid- vapor  meniscus.  For 
the  wickless  micro  heat  pipe,  ATcr;t  is  the  critical  temperature  drop  across  the  liquid 
film.  Since  Ts  =  Tv  +  Arcrit,  eqn.  (1.7)  can  be  written  as 
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The  maximum  heat  flux  related  to  the  boiling  limit,  assuming  heat  conduc¬ 
tion  is  dominant  in  the  thin  liquid  film,  is  given  by 


<7crit  —  ATcrit/  & 


(1.9) 
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It  can  be  seen  from  eqns.  (1.8)  and  (1.9)  that  the  critical  heat  flux  strongly 
depends  on  the  heat  pipe  working  condition.  The  dominant  term  in  eqn.  (1.8)  is 
2cr ( 1  / r*  —  1  /Rm)/pv.  For  a  fixed  r^,  a  higher  vapor  temperature  results  in  a  much 
higher  vapor  pressure  pv,  and  therefore  a  much  lower  critical  heat  flux.  Since  for  many 
mainframe  computers,  the  maximum  operating  junction  temperature  for  functionality 
is  about  80°  C,  the  boiling  limit  was  examined  under  this  working  temperature.  Figure 
1.6  shows  the  boiling  limits  as  a  function  of  liquid  thickness  6  at  T^  =  80° C  and  r j,  = 
2.54  x  10~7m.  For  water  as  the  working  fluid,  the  critical  heat  flux  is  generally  high 
for  most  of  the  range  of  6  in  the  figure. 

For  a  larger  liquid  film  thickness  (T>  5  x  10~5  m),  however,  the  critical  heat  flux 
drops  below  100  W/cm2.  For  methanol  as  the  working  fluid,  the  critical  heat  flux 
is  generally  low.  When  the  film  thickness  becomes  larger,  the  critical  heat  flux  even 
drops  to  below  10  W/cm2.  In  the  above  calculation,  the  radius  of  the  liquid- vapor 
meniscus  was  assumed  to  be  much  larger  than  r^,  and  its  effect  on  the  critical  heat 
flux  was  neglected.  For  high  temperature  heat  pipes  with  liquid  metal  working  fluids, 
the  boiling  limit  rarely  occurs  (Cao  and  Faghri,  1992b),  so  it  was  not  considered  here. 

The  boiling  limit  also  largely  depends  on  the  effective  radius  of  small  vapor 
or  noncondensible  gas  bubbles,  rj.  Figure  1.7  presents  the  boiling  limits  as  a  function 
of  rf,  for  Tv  =  80°  C  and  6  =  10-5  m.  Both  critical  heat  fluxes  for  water  and  methanol 
drop  to  low  values  for  relatively  larger  r&,  even  for  a  relatively  thin  liquid 
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Figure  1.6:  Boiling  limits  as  a  function  of  liquid  thickness  for  Tv  =  80°C,  h 
2.54  x  10-7  m  (Cao  and  Faghri,  1994) 


film  thickness.  To  reduce  the  value  of  rj,  complete  cleaning  and  degassing  of  the  heat 
pipe  inner  wall  surface  are  needed.  Since  the  vapor  space  of  a  micro  heat  pipe  is  very 
small,  difficulties  related  to  the  heat  pipe  cleaning  and  degassing  are  expected. 


1.4  Capillary  Limitation 

Like  conventional  heat  pipes,  the  capillary  limit  remains  one  of  the  major  operating 
limits  for  miniature  and  micro  heat  pipe  designs.  Cotter  (1984)  pointed  out  that  the 
heat  transport  capacity  of  micro  heat  pipes  depends  more  critically  on  the  volume  of 
working  fluid  available  than  conventional  heat  pipes  with  wicks.  The  mean  radius  of 
curvature  R  is  given  by 


R  R] 

where  Ri  and  f?2  are  two  principal  radii  of  curvature  for  the  surface.  If  the  surface 
tension  is  denoted  by  <7,  and  the  pressures  in  the  liquid  and  vapor  are  pt  and  then 


Pv~Pt  = 


R 


(1.11) 


d  .  .  d  f  a 

~  p,)  =  TAr 


(1,12) 


The  approximate  maximum  heat  transport  for  steady-state  derived  by  Cot- 


21 


ter  (1984),  including  the  use  of  eqn.  (1.12),  is 


Qmax  — 


0.1  6P(K£/Kv)1/2\  (ahjg\  /^y/2  / A0-75 
8  7 r  H(  1)  I  V  vt  )  \uv)  l  Lt 


(1.13) 


where  IQ  and  I\v  are  the  flow  shape  factors,  H(l)  is  an  integral  over  the  length  of 
the  pipe  as  given  by  Cotter  (1984),  A  is  the  total  cross  sectional  area  of  the  heat 
pipe,  and  0  is  a  dimensionless  geometrical  factor.  As  an  illustrative  case,  Cotter 
(1984)  considered  a  micro  heat  pipe  with  an  equilateral  triangular  cross  section  filled 


with  the  optimum  amount  of  methanol  to  operate  at  50°C  with  a  uniformly  heated 


evaporator  and  uniformly  cooled  condenser.  For  a  triangular  pipe  with  a  side  of  0.02 
cm  and  a  length  of  1  cm,  the  maximum  heat  transport  Qma.x  was  found  to  be  0.03  W. 
If  placed  closely  together  on  a  surface,  micro  heat  pipes  of  this  size  could  provide  a 
few  watts  per  square  centimeter  of  cooling.  Cotter  concluded  that  micro  heat  pipes 
are  candidates  for  thermal  control  of  micromechanical  and  microelectronic  devices  in 
cases  where  thermal  loads  are  modest,  and  where  a  high  standard  of  temperature 
constancy  and  uniformity  are  required. 

Babin  et  al.  (1989)  compared  the  maximum  heat  transport  (Qmax,  eqn. 
(1.13))  from  Cotter  (1984),  with  their  experimental  data  related  to  the  capillary  limit. 
They  found  that  eqn.  (1.13)  underpredicted  the  heat  transport  limit  substantially, 


especially  when  the  operating  temperature  was  relatively  high.  The  maximum  heat 


transport  relation  based  on  the  conventional  larger  heat  pipe  was  also  presented 
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and  compared  with  related  experimental  data  by  Babin  et  al.  (1989).  Unlike  the 
micro  heat  pipe,  the  conventional  method  involves  defining  vapor  pressure  and  liquid 
pressure  separately,  and  then  combining  them  to  obtain  the  maximum  heat  flux.  The 
relation  obtained  by  Babin  et  al.  (1989)  based  on  Chi  (1976)  is 

2<r/reff  -  ptgLt  sin  ip 
Pi  C  (fvRev)  pv 

K  AghjgPi  2(rh,v)2Avpvhfg 

where  <p  is  the  heat  pipe  tilt  angle,  K  is  the  wick  permeability,  rh^v  is  the  hydraulic 
radius  of  the  vapor  space,  and  C  is  a  modification  coefficient.  The  values  of  C  and 
fvRev  depend  on  the  vapor  flow  regime.  Good  agreement  was  obtained  between  eqn. 

(1.14)  and  the  related  experimental  data  (Babin  et  ah,  1989).  The  agreement  may 
be  due  to  the  relatively  large  size  of  the  heat  pipe  tested.  Experimental  data  related 
to  actual  micro  heat  pipes  are  needed  for  further  examination  of  eqns.  (1-13)  and 

(1.14) . 

Gerner  et  al.  (1992)  derived  a  correlation  for  the  maximum  heat  transport 
in  the  form 


max  — 

Teff 


(1.14) 


n  —  n  ^  a  fill  r)3 

V  max  —  ^  r-  ~  ^  -Ls  h 


(1,15) 


2048  vvLt 

The  constant  C  was  determined,  using  the  experimental  data  from  Babin  et  al.  (1989), 
to  be  2.38.  Therefore,  eqn.  (1.15)  can  be  written  as 


23 


ahfgDl 


(1.16) 


Qmax  =  0.01 

vvLt 

For  a  water  heat  pipe  with  hydraulic  diameter  of  100  pm  and  length  of  2.54  cm,  for 
example,  heat  transport  on  the  order  of  1  mW  can  be  expected  per  device.  If  127 
of  these  micro  heat  pipes  per  inch  are  used,  and  assuming  that  all  of  the  energy  is 
transferred  through  the  end  of  the  evaporator,  a  heat  flux  up  to  15  W/cm2  can  be 
dissipated  under  a  working  temperature  of  70°C  (Gerner  et  ah,  1992). 

Since  micro  heat  pipes  are  very  small,  the  evaporation  in  the  liquid  “micro¬ 
layer”  region  may  play  a  role  during  heat  transfer  in  micro  heat  pipes.  The  microlayer 
region  is  also  referred  to  as  the  contact  line  region.  A  thin  film  of  liquid  on  a  horizontal 
solid  surface  is  schematically  shown  in  Fig.  1.8.  For  a  very  thin  film  on  the  surface, 
the  pressure  gradient  for  the  liquid  flow  is  due  to  the  conventional  capillary  force  and 
the  disjoining  pressure,  which  results  from  the  repulsion  of  the  vapor  phase  by  the 
solid  and  the  liquid  from  long  range  intermolecular  forces.  The  chemical  potential 
of  an  evaporating  thin  film  differs  from  that  of  the  bulk  liquid  by  an  amount  often 
referred  to  as  the  excess  potential  due  to  the  solid/liquid/vapor  intermolecular  force 
field.  This  potential  is  defined  in  units  of  pressure  (disjoining  pressure).  The  pressure 
change  at  the  liquid-vapor  interface  is  (Schonberg  and  Wayner,  1990;  Dasgupta  et 
ah,  1991) 
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Figure  1.8:  Evaporating,  nonevaporating  regions  and  coordinate  system  (not  to  scale) 
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(1.17) 


Pv  -  Pi  -  p  +  crKc 

where  A  =  A/6tt ,  A  is  the  Hamaker  constant,  a  is  the  surface  tension,  and  Kc  is  the 
curvature  of  the  interface.  When  the  liquid  film  is  thin  with  a  relatively  flat  surface, 
the  second  term  in  eqn.  (1.17)  can  always  be  neglected  compared  to  the  first  term.  For 
a  very  thin  film,  8  <  8 o,  where  So  is  the  nonevaporating  film  thickness  shown  in  Fig. 
1.8,  a  thin  film  region  persists  at  equilibrium  with  its  saturated  vapor,  even  though 
it  is  superheated.  The  film  does  not  evaporate  due  to  the  reduced  vapor  pressure  of 
the  film  caused  by  the  disjoining  pressure.  This  region  is  called  the  nonevaporating 
region  where  8  <  80.  Beyond  this  region  is  the  microlayer  region  where  evaporation 
and  liquid  flow  occur.  The  mass  flux  of  vapor  leaving  the  liquid/ vapor  interface  is 
modeled  as  (Schonberg  and  Wayner,  1990) 

m' =  a  (T£v  -  Tv)  +  b  (pi  -  pv)  (1-18) 

where 


a  =  2.0 


/  M  \1/2 

\2ttRT£v) 


( pvM  hfg\ 

\  rtvt£v  J 


6  =  2.0 


M 

2-7T  RT£v 


1/2 


ViPv  \ 
RT£vJ 
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T(v  is  the  temperature  of  the  liquid  at  the  liquid-vapor  interface,  and  Vi  is  the  molar 
volume  of  the  liquid.  A  closed-form  solution  for  the  dimensionless  integral  contact 
line  heat  sink  was  obtained  by  neglecting  the  interface  curvature  and  applying  the 
lubrication  theory  related  to  the  mass  flow  in  the  thin  film  (Schonberg  and  Wayner, 
1990) 


Q*  =  [18(1  +  A3)ln[^1+.A)]  +  6(tT3  -  1)  +  9A(1  -  r2)  +  18  A -  l)]1/2  (1.19) 

l  +  AT] 

The  integral  heat  sink  for  the  contact  line  region,  which  is  defined  as  the  heat  absorbed 
by  the  microlayer  region  between  r)  =  1  and  rj,  is  directly  related  to  Q*,  which  is 

Q^Ahl^Q-  (1.20) 

oIS^Xq 

where 


x  =  (Lhih  (1.21) 

k 

<1=J  (1-22) 

<20 


A 

VidAT 


(1.23) 
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bA' 

aAT 


1/3 


(1.24) 


*o  = 


A T  =  TS-  Tv  (1.25) 

Figure  1.9  presents  Q*  as  a  function  of  rj  for  sodium,  water  and  methanol  as  heat 
pipe  working  fluids.  For  relatively  large  values  of  rj  =  6 /So,  Q*  approaches  a  constant, 
which  is  the  maximum  dimensionless  contact  line  heat  sink  Qmax-  Qmax  corresponds 
to  a  maximum  dimensionless  liquid  film  thickness,  r/max,  in  the  microlayer  region  in 
which  the  capillary  pressure  is  negligible.  The  corresponding  maximum  liquid  film 
thickness  is 


^max  —  7/max^O  (1.26) 

From  Fig.  1.9,  the  maximum  thickness  of  the  microlayer  is  on  the  order  of  0.5  ~  1.0 
fim.  For  a  micro  heat  pipe  with  hydraulic  radius  on  the  order  of  10  \i m,  the  liquid 
film  thickness  in  the  axial  direction  close  to  the  evaporator  end  cap  should  be  on  the 
order  of  1  fi m,  which  is  in  the  microlayer  region.  Therefore,  the  disjoining  pressure 
may  help  to  draw  liquid  toward  the  evaporator  end  cap  and  increase  heat  transport 
capacity.  In  the  circumferential  direction,  as  shown  in  Fig.  1.10  for  a  micro  heat  pipe 
with  a  triangular  cross  section,  the  disjoining  pressure  may  help  to  spread  liquid  from 
the  corners  to  the  rest  of  the  inner  wall  surface  and  improve  the  cooling  condition  at 
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the  evaporator.  As  a  result,  the  heat  transport  capacity  for  a  micro  heat  pipe  may  be 
somewhat  larger  than  the  capillary  limit  derived  based  on  conventional  larger  heat 
pipes.  It  should  be  pointed  out  that  the  above  arguments  are  based  on  a  simplified 
theory,  and  are  qualitative  in  nature.  A  detailed  experimental  study  is  needed  for 
estimating  the  increase  in  the  heat  transport  capacity  quantitatively. 

1.5  Conclusions  and  Remarks 

The  literature  review  revealed  that  the  miniature  flat  heat  pipes  developed  by  Plesch 
et  al.  (1991)  and  other  investigators  were  capable  of  removing  heat  fluxes  as  high 
as  60  W/cm2  under  the  thermal  siphon  operating  mode.  The  maximum  heat  flux 
obtainable  may  be  further  increased  by  optimizing  the  longitudinal  grooves  on  the 
inner  flat  surface.  Further  performance  improvements  can  also  be  made  by  reducing 
the  contact  resistance  between  the  heat  pipe  and  the  electronic  chip.  Since  micro 
electronic  components  are  shaped  as  thin  plates  with  flat  base  surfaces,  and  the  heat 
resistance  in  the  transverse  direction  is  small  due  to  the  high  thermal  conductivity  of 
silicon  (on  the  order  of  100  W/(m-K)),  most  of  the  heat  generated  in  the  electronic 
component  can  be  transferred  to  the  flat  heat  pipe  and  dissipated  via  the  heat  pipe 
condenser  without  causing  large  temperature  gradients  in  the  electronic  component. 
For  the  current  highest  chip  heat  fluxes  on  the  order  of  60  W/cm2,  these  miniature 
flat  heat  pipes  can  be  expected  to  have  a  large  commercial  application. 
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As  the  size  of  the  heat  pipe  decreases,  the  micro  heat  pipe  may  encounter 
the  vapor  continuum  limitation.  This  limitation  may  prevent  the  micro  heat  pipe 
from  working  under  lower  temperatures.  For  a  heat  pipe  operating  at  higher  working 
temperatures,  this  limitation  is  temporary  and  will  disappear  when  the  heat  pipe 
temperature  increases.  In  addition  to  the  vapor  continuum  limitation,  the  micro  heat 
pipe  is  also  subject  to  the  operating  limits  of  the  conventional  heat  pipe.  Of  these 
operating  limits,  the  capillary  limit  remains  the  most  important  for  the  micro  heat 
pipe.  The  maximum  heat  flux  based  on  the  conventional  formulation  is  on  the  order 
of  10  W/cm2,  which  is  relatively  low.  However,  since  the  size  of  micro  heat  pipes  is 
very  small,  the  disjoining  pressure  may  help  to  increase  the  heat  transport  capacity 
of  the  heat  pipe.  Further  study  is  needed  to  quantify  this  effect.  Finally,  in  order  for 
micro  heat  pipes  to  find  commercial  application  in  microelectronic  cooling,  it  must 
compete  with  other  cooling  methods  such  as  forced  convection,  impingement  and 
two-phase  direct  cooling  in  areas  such  as  manufacturing  cost  and  reliability.  Because 
of  the  extremely  small  dimension  of  micro  heat  pipes,  difficulties  in  manufacturing 
and  the  subsequent  cleaning  and  degassing  are  expected.  Furthermore,  the  “true” 
micro  heat  pipe  built  directly  into  the  silicon  must  compete  with  the  miniature  heat 
pipe  attached  on  the  chip  surface  in  this  area,  which  is  much  less  costly.  Also,  the 
operating  conditions  of  the  miniature  heat  pipe  could  be  further  improved  by  reducing 
the  contact  resistance  between  the  heat  pipe  and  the  electronic  chip. 
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Chapter  2 


THERMAL  ANALYSIS  OF  A 
MICRO  HEAT  PIPE 


2.1  Summary 

A  detailed  mathematical  model  is  developed  in  which  the  heat  and  mass  transfer  pro¬ 
cesses  in  a  micro  heat  pipe  (MHP)  are  examined.  The  model  describes  the  distribution 
of  the  liquid  in  an  MHP  and  its  thermal  characteristics  depending  upon  the  liquid 
charge  and  the  applied  heat  load.  The  liquid  flow  in  the  triangular-shaped  corners  of 
an  MHP  with  polygonal  cross  section  is  considered  by  accounting  for  the  variation  of 
the  curvature  of  the  free  liquid  surface  and  the  interfacial  shear  stresses  due  to  liquid- 
vapor  frictional  interaction.  The  predicted  results  obtained  are  compared  to  existing 
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experimental  data.  The  importance  of  the  liquid  fill,  minimum  wetting  contact  angle, 
and  the  shear  stresses  at  the  liquid- vapor  interface  in  predicting  the  maximum  heat 
transfer  capacity  and  thermal  resistance  of  the  MHP  is  demonstrated. 

2.2  Introduction 

In  an  MHP,  condensate  return  is  accomplished  by  the  flow  of  liquid  in  the  corners 
of  the  polygonal  heat  pipe  container;  the  condensate  being  held  in  the  corners  by 
surface  tension.  The  objective  of  the  present  study  is  to  develop  a  mathematical 
model  which  can  determine  the  maximum  heat  capacity  and  thermal  resistance  of 
MHPs.  The  present  model  incorporates  the  following  new  features: 

•  The  influence  of  the  liquid  charge  on  the  length  of  the  liquid  blocking  zone  and 
on  the  liquid  distribution. 

•  The  heat  transfer  through  a  liquid  film  in  the  evaporator  is  described  with 
respect  to  the  disjoining  pressure. 

•  The  heat  transfer  through  a  film  of  condensate  in  the  condenser  is  considered 
and  the  closed-form  solutions  are  obtained  using  an  approach  similar  to  that 
developed  for  axially-grooved  heat  pipes. 

•  The  influence  of  shear  stresses  at  the  free  liquid  surface  in  a  corner  due  to  the 
frictional  vapor-liquid  interaction  on  the  liquid  flow  is  taken  into  consideration. 
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•  The  possibility  of  the  variation  of  the  meniscus  contact  angle  in  the  condenser 
is  treated  by  considering  the  condensate  film  formation. 

•  The  model  describes  MHP  operation  in  the  case  of  convective  heat  transfer 
outside  the  evaporator  and  condenser  as  well  as  that  for  a  prescribed  heat  flux 
distribution. 


2.3  Mathematical  Model 

The  present  model  deals  with  a  polygonal  MHP  having  N  triangular  corners,  which 
can  easily  be  extended  for  cross  sections  having  different  corner  configurations.  Figure 
2.1  presents  cross  sections  of  the  triangular  micro  heat  pipe  used  to  illustrate  the 
methodology  of  the  present  study.  The  heat  load  is  uniformly  distributed  between 
all  of  the  corners.  Near  the  evaporator  end  cap,  if  the  heat  load  is  sufficient,  the 
liquid  meniscus  is  depressed  in  the  corner,  and  its  cross  sectional  area  and  radius 
of  curvature  of  the  free  surface  are  extremely  small.  Most  of  the  wall  is  dry  or  is 
covered  by  a  nonevaporating  liquid  film.  In  the  adiabatic  section,  the  liquid  cross 
sectional  area  is  comparatively  larger.  The  inner  wall  surface  may  be  covered  with  a 
thin  liquid  film  due  to  the  disjoining  pressure.  At  the  beginning  of  the  condenser,  a 
film  of  condensate  is  present  on  the  wall,  and  liquid  flows  through  this  film  toward 
the  meniscus  region  under  the  influence  of  surface  tension.  For  normal  operation,  the 
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contact  angle  of  the  meniscus  6  is  constant  and  equal  to  the  minimum  wetting  contact 
angle  6q  in  the  evaporator  and  adiabatic  sections  and  the  beginning  of  the  condenser 
zone,  which  is  fixed  for  a  specific  working-fluid/ container  combination.  In  Fig.  2.1, 
the  possibility  of  liquid  blocking  the  end  of  the  condenser  is  shown.  The  details  of 
the  evaporator  and  condenser  sections  are  shown  in  Fig.  2.2. 

The  fluid  flow  and  heat  transfer  within  a  micro  heat  pipe  operating  under 
steady-state  conditions  are  modeled,  neglecting  axial  heat  conduction  in  the  wall  and 
liquid.  The  conservation  equations  for  the  steady-state  operation  of  an  MHP  are  the 
continuity,  momentum  and  energy  equations  for  the  liquid,  vapor  and  wall,  and  the 
Laplace- Young  equation  for  the  radius  of  curvature  at  the  liquid- vapor  interface.  At 
any  axial  location,  the  following  mass  conservation  equation  must  hold  over  the  cross 
section  of  the  MHP: 

wvpvAv  =  NwipiAi  =  Nriie  (2.1) 

where  rri(  is  the  mass  flow  rate  through  a  corner.  The  average  liquid  and  vapor 
velocities  in  the  axial  direction  in  Eq.  2.1  are: 

Wi=^-[  [  we(x,y)  dxdy,  wv  =  —  [  [  wv{x,  y)  dx  dy  (2.2) 

A-l  J  J Ai  -Ai-v  J  J Av 

A  MHP  contains  a  definite  amount  of  a  working  fluid  Mt,  which  is  dis¬ 
tributed  in  accordance  to  the  following  relation: 
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Figure  2.2:  Details  of  the  micro  heat  pipe  cross  section  in  (a)  the  evaporator,  (b)  the 
condenser 
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rLt—Lb  i-Lt-Lt, 

Mt  =  N  piAi  dz  + 

Jo  Jo 


pvAv  dz  +  Mb  +  M$ 


(2.3) 


where  Mh  and  Lb  are  the  mass  and  length  of  the  liquid  which  blocks  the  condenser, 
and  Ms  is  the  mass  of  the  thin  liquid  films  in  the  condenser  and  adiabatic  sections  in 
the  interval  0  <  s  <  Lx  (see  Fig.  2.2(b))  for  all  of  the  MHP  walls,  while  At  doesn’t 
include  the  thin  film  existing  in  this  interval.  It  is  anticipated  that  the  excessive 
liquid,  the  amount  of  which  depends  upon  the  operating  conditions,  is  contained  at 
the  end  of  the  condenser  because  of  the  MHP  small  inner  diameter  and  the  effect  of 
the  vapor  flow  on  the  liquid  (Fig.  2.1). 

The  axial  conservation  of  momentum  equation  for  an  incompressible  vapor 
flow  using  the  one-dimensional  boundary-layer  approximation  is: 

-j ~(Py  +  pvgz  sirup  +  pvfivwl)  =  ~  (2-4) 

UZ  -*-dh,v 

where  8V  is  the  momentum  flux  coefficient,  fv  is  the  friction  coefficient,  and  Dh,v  is 
the  hydraulic  diameter  of  the  vapor  channel.  The  values  of  these  coefficients  can 
be  determined  using  the  results  of  previous  investigators  for  the  two-dimensional, 
laminar  boundary-layer  equations  with  suction  and  injection.  For  example,  for  the 
small  radial  Reynolds  numbers  usually  seen  in  MHPs  (0  <  Rer  <  2),  in  the  evapo¬ 
rator  and  adiabatic  sections,  the  results  given  by  Bankston  and  Smith  (1972)  can  be 
approximated  as  follows: 
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Pv  =  1.33  —  0.0075Rej? 


x  fl6  +  0.25Re^ 

/v  V  Re  J’ 

In  the  condenser  section  (Bowman  and  Hitchcock,  1988;  Jang  et  ah,  1991): 


fv  =  16[1.2337  -  0.2337exp(— 0.0363Rer)][exp(1.2Ma)]/Re,  f3v  =  1.33 

The  velocity  of  the  liquid  phase  is  very  small  in  comparison  to  that  of  the  vapor 
flow,  so  the  interfacial  shear  for  the  vapor  is  computed  by  assuming  the  liquid  to  be 
stationary  (Longtin  et  ah,  1992). 

The  axial  transport  of  condensate  in  an  MHP  takes  place  in  the  corners, 
where  most  of  the  liquid  resides.  The  thickness  of  the  liquid  films  in  the  condenser 
and  adiabatic  sections  in  the  interval  0  <  s  <  L\  (see  Fig.  2.2(b))  is  assumed  to  be 
too  small  to  contribute  to  the  axial  mass  transport.  Since  the  axial  Reynolds  number 
for  liquid  flow  in  an  MHP  is  usually  less  than  1,  it  can  be  considered  to  be  viscous 
and  quasi-one-dimensional.  Thus,  the  conservation  of  momentum  equation  for  the 
liquid  flow  in  a  corner  with  cross  sectional  area  At  is: 


dP, 


,  Pi  d  2  a  \  ,  •  (  d2we  ,  d2we\ 

+  —  +  pi9  sm  r  =  ^  ~jgr  +  -Q^) 


(2.5) 


dz  2  At  dz 

The  second  term  in  Eq.  2.5  is  the  contributions  of  the  dynamic  component  of  the  flow 
rate  and  the  axial  variation  of  the  cross  sectional  area  on  the  liquid  pressure  gradient. 
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The  present  numerical  results  have  shown  that  this  term  changes  the  liquid  pressure 
gradient  by  less  than  0.1  percent. 

The  conservation  of  energy  equation  is  cast  into  a  form  which  reflects  the 
change  in  the  axial  mass  flow  rate  of  liquid  due  to  evaporation  and  condensation: 

+  \  o <z<Le 

hfg  V  tle  / 

0  (adiabatic  section)  Le  <  z  <  Le  +  La  (2-6) 

7 — heXiC(Tex,c  —  Tv)  (l  +  -=^')  ,  Le  +  La  <  z  <  Lt  —  Lb 

hfg  \  hc  J 

where  he  and  hc  are  the  mean  internal  heat  transfer  coefficients  between  the  wall 
and  vapor,  and  hex,e  and  hex,c  are  the  heat  transfer  coefficients  between  the  wall  and 
the  ambient.  The  mean  coefficients  he  and  hc  for  evaporation  and  condensation, 
respectively,  include  the  resistances  of  the  wall  and  liquid  film,  which  are  given  by  the 
heat  transfer  analysis  for  each  of  the  heat  loaded  sections  at  every  point  on  z.  The 
vapor  temperature  is  denoted  by  Tv.  For  a  given  axial  heat  load  function  Q(z),  Eq. 
2.6  can  be  rewritten  as: 


—(ptW(Ae) 

dz 


1  dQ(z) 
h  fgN  dz 


(2.7) 


For  the  case  of  radiative  heat  transfer  in  the  condenser  region,  the  heat  transfer 


coefficient  in  Eq.  2.6  is: 
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hex,c  —  coc(Tv  +  TeXtC)(T^  +  TlJ  (2-8) 

An  iterative  procedure  can  be  used  to  determine  hex>c  from  Eq.  2.8,  since  it  is  a  strong 


function  of  Tv  in  the  case  of  a  radiative  boundary  condition. 

The  interfacial  meniscus  radius  of  curvature  is  related  to  the  pressure  dif¬ 
ference  between  the  liquid  and  vapor  by  the  Laplace- Young  equation,  which,  in  dif¬ 
ferential  form,  is: 


dPi  dPv  d  f  a  \ 
dz  dz  dz\Rm ) 

The  vapor  temperature  and  pressure  are  related  by  the  perfect  gas  equation: 


(2.9) 


Pv  =  PvRgTy  (2.10) 

Equation  2.10  is  accurate  enough  for  comparatively  low  pressure  values  for  the  condi¬ 
tions  considered  in  the  present  study  (the  maximum  error  value  was  0.73%).  However, 
in  the  general  case,  one  should  use  an  empirical  equation  based  on  saturation  condi¬ 
tions  and  the  Clausius- Clapeyron  equation  instead  of  Eq.  2.10. 

The  mean  axial  liquid  and  vapor  velocities  at  the  evaporator  and  condenser 
end  caps  are: 


wi\z=0  =  Wv\z= 0  =  0 


(2.11) 
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Referring  to  Fig.  2.1,  the  liquid  in  the  corner  is  attached  to  two  walls,  and 
interacts  with  the  vapor  flow  at  the  liquid- vapor  interface,  u>2-  Therefore,  the  no-slip 
condition  and  relation  for  the  equality  of  tangential  shear  stress  on  these  surfaces  are: 


=  0; 


dwt 


dn 


U>2 


fv  — 2 

'2ntPvWv 


>.13) 


The  second  boundary  condition  in  Eq.  2.13  is  an  essential  feature  of  the  present 
model.  Here  n  is  the  normal  vector  on  the  free  surface  of  the  liquid  meniscus.  The 
cross  section  of  the  vapor  channel  has  a  complicated  shape  and  the  assumption  of  a 


circular  cross  section  is  not  completely  exact  for  N  <3,  but  it  is  very  close  for  N  >  3, 
where  the  vapor  channel  is  nearly  circular. 

The  vapor  and  liquid  pressures  at  the  evaporator  end  cap  are: 


Py\z=Q  —  PvOi  Pt\z=0  —  PvO - p  (2-14) 

The  values  of  Pv 0  and  Rm 0  are  to  be  determined  using  additional  conditions.  When 
solving  Eq.  2.6,  where  Tv(z)  is  also  unknown  and  depends  on  Pv,  the  value  of  Pv0  must 
be  chosen  such  as  to  satisfy  Eq.  2.12.  When  Q(z)  and  Tv0  are  known,  Pv 0  =  Pv(Tv0), 
which  can  be  obtained  using  the  perfect  gas  law,  Eq.  2.10.  The  radius  of  curvature  of 
the  meniscus  at  the  evaporator  end  cap,  Rm 0,  is  found  using  Eq.  2.3,  since  the  value 


of  Rm o  influences  the  axial  distribution  of  liquid.  In  general,  the  radius  of  curvature 
of  the  meniscus  is  bounded  by  i?m, min  <  Rm o  <  Rm, max-  When  liquid  blocks  part  of 
the  condenser  section,  Rm0  is  specified  such  that  Rm  ~  i?m,max  at  the  location  where 
the  liquid  blocking  begins,  where  i?m,max  is  approximately  the  radius  of  the  largest 
circle  which  can  be  inserted  in  the  inner  MHP  cross  section,  (Fig.  2.1).  For  the  case 
in  which  no  liquid  blocks  the  condenser,  the  value  of  Rm0  is  given  by  Eq.  2.3  where 
—  0  and  —  0,  or  RmO  —  Rm,mirL  when  CJa  —  Qmax- 

To  determine  the  heat  transfer  coefficients  ht  and  hc  in  Eq.  2.6,  the  heat 
transfer  during  evaporation  and  condensation  in  a  triangular  capillary  groove  is  ex¬ 
amined. 


2.4  Heat  Transfer  in  the  Condenser  Section 


The  equation  for  the  condensate  film  thickness  distribution,  6(s),  shown  in  Figs.  2.1 
and  2.2,  is  given  by  Kamotani  (1976a)  for  a  low-temperature  heat  pipe  with  axial 
grooves: 


PeaP  =  h{Tv  -  Tw)  r>  ds 
ds3  hjg  Jo  6 

For  the  case  of  a  polygonal  MHP,  the  boundary  conditions  are: 


(2.15) 
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d6  _Q  £1  -0-  —  -  — • 

ds  s= 0  ’  ds 3  s=0  ’  ds2  s-Lx  Rm' 


ds 


=  tan  9 

s—L\ 


(2.16) 


Equation  2.15  has  been  derived  from  the  conservation  of  mass,  momentum  and  energy 
equations  for  the  flow  of  a  thin  film  and  the  relation  between  the  thickness  and  the 
curvature.  The  interfacial  thermal  resistance  is  usually  neglected  in  comparison  to 
the  resistance  of  the  liquid  film,  but  it  can  be  accounted  for  in  the  expression  for 
the  condenser  heat  resistance.  Equation  2.15  can  be  solved  for  the  film  thickness 
by  numerical  methods.  In  the  present  paper,  however,  an  approximate  closed-form 
solution  is  obtained  using  the  assumption  that  the  liquid  film  thickness  variation  along 
the  s-direction  is  small  in  comparison  to  its  thickness;  i.e.,  the  film  is  assumed  to  be 
nearly  planar.  The  film  thickness  in  the  middle  region  of  the  wall  (see  Fig.  2.2(b)) 
can  be  approximated  by: 


8_ 

B 


=  Co  +  Cl—  +  c2—  +  c3—  +  c4— 


From  the  boundary  conditions  in  Eq.  2.16,  the  coefficients  are: 


(2.17) 


Ci  =  C3  =  0;  C4  = 


Li/Rm  —  tan# 
8(Ti/ B)z  ’ 


r  b  rxn2 

Ci  -  - 0  ”77  C4 


The  value  of  Co  is  obtained  by  differentiating  Eq.  2.15,  inserting  Eq.  2.17  and 


leglecting  the  term  which  contains  the  first  derivative  of  8.  The  mean  heat  transfer 
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coefficient  of  the  thin  film  region,  hcj  =  £p/<5,  is  approximated  as: 


lc,f 


BC0 


hfgpecrkffa/Rm  -  tan  9) 


1  0.25 


tHL\(Tv  -  Tw) 


(2,19) 


From  Fig.  2.2(b),  the  value  of  Li  is  defined  as: 


£l  =  ^  -  -  7-^—B 

2  sin  7  tan  7 


(2.20) 


Equations  2.19  and  2.20  are  solved  for  Co  and  Lx  using  an  iterative  procedure.  The 
increase  of  the  liquid  film  thickness  along  the  s  coordinate  in  the  interval  5  =  0  to 
5  =  Lx  in  the  numerical  experiments  is  less  than  8.0  percent,  so  Eq.  2.19,  where  it  is 
assumed  that  6  =  BCo ,  gives  hcj  to  within  4.0  percent. 

The  heat  transfer  coefficient  in  the  meniscus  region  can  be  estimated  as 
follows.  For  s  >  Li,  the  film  thickness  increases  sharply  (see  Fig.  2.2(b)): 


&  —  Bm  + 


R2m  +  (s  —  Lj)2  +  2(5  —  Li)Rm  sin  9 


0.5 


(2.2i; 


Then,  for  the  mean  heat  transfer  coefficient  in  the  meniscus  region,  hc,men,  we  have: 


-L 


B /2  kt  ds 


c,men  )Ll  (B/2  -Lx)8 
For  6  >  20°  there  is  a  simple  analytical  approximation  for  hC)men: 


(2.22) 
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6  =  6\s=Ll+sf  (2.24) 

as  s=h 

The  mean  heat  transfer  coefficient  in  the  condenser  section  for  a  given  z  can 
be  obtained  by  combining  hcj  and  hc<men  with  their  respective  areas,  and  including 
the  thermal  resistance  of  the  wall: 


K =  ([ir*'"'  +  l1  -  ]_1+|^)_1  (2-25) 

The  wall  temperature  in  Eq.  2.19  is  determined  by  the  thermal  resistance  between 
the  vapor  temperature  and  the  ambient  temperature: 

Tw  =  Ter,c  +  (Tv  -  TeXiC)^—(  -±-  +  iV1  (2.26) 

Therefore,  Eqs.  2.19  -  2.26  are  solved  for  hc  and  Tw  using  an  iterative  procedure. 

Depending  on  the  amount  of  working  fluid  within  the  MHP,  part  of  the 
condenser  may  be  blocked  with  liquid.  In  this  case,  the  effective  condenser  length  is 
Lc  —  Lb,  where  Lb  is  the  length  blocked  by  liquid. 
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(2.27) 


Lb  = 


Mh 


Mb 


(Av  +  N  At)  pi  At  pi 

The  overall  thermal  resistance  of  the  condenser  can  then  be  determined  as  follows: 


rr  = 


i _ / 

'  r  Lh  1  J  ( . 


(M>)  dQ  11, 
—  dz 


(2.28) 


Qa(dJC  Lb)  J (Le+La)  dz  BA  hc 
In  the  numerical  procedure,  for  the  case  of  convective  heat  transfer  outside  the  MHP, 

dQ / dz  is  defined  using  Eq.  2.7. 

Now,  the  variation  of  the  contact  angle  of  the  meniscus  along  the  condenser 
is  to  be  considered.  At  the  inlet  of  the  condenser,  the  contact  angle  is  6  —  90.  In 
the  region  near  the  liquid  blocking  zone,  9  =  0.  Therefore,  variation  of  the  contact 
angle  within  the  condenser  takes  place.  At  the  beginning  of  the  condenser,  usually 
Li/Rm  >  tan#  and  the  liquid  film  is  very  thin.  Along  the  condenser,  while  Li/Rm  is 
decreasing,  there  is  a  z  location  where  L\/ Rm  ~  tan  9.  In  the  vicinity  of  this  location, 
Eq.  2.15  cannot  be  solved;  that  is,  the  conservation  of  mass  and  energy  in  the  film 
cannot  be  satisfied.  This  causes  the  contact  angle  9  to  begin  changing,  starting  from 
this  particular  z  location  in  the  condenser.  In  the  present  model,  it  is  assumed  that 
the  9  variation  after  this  location  along  the  z-axis  in  the  condenser  is  described  by 
the  relation: 


Rm 


tan  9  =  0.01 


(2.29) 
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The  value  0.01  is  taken  from  both  physical  and  numerical  considerations.  The  choice 
of  this  value  is  made  since  it  is  very  close  to  zero  in  comparison  to  that  given  by  the 
right  hand  side  of  Eq.  (2.29)  at  the  beginning  of  the  condenser,  which  can  be  greater 
than  unity.  This  value,  unlike  zero,  still  provides  a  physically  reasonable  solution  for 
Eqs.  2.15  and  2.16.  This  method  allows  for  the  approximate  equality  of  the  contact 
angles  of  the  liquid  flow  in  the  corner  and  at  the  beginning  of  the  liquid  blocking 
zone. 


2.5  Heat  Transfer  in  the  Evaporator  Section 

Following  the  approaches  of  Stephan  and  Busse  (1992),  Solov’yev  and  Kovalev  (1984) 
and  Holm  and  Goplen  (1979)  for  the  description  of  the  heat  transfer  in  the  meniscus/ thin- 
film  transition  region,  an  approximate  thermal  analysis  of  evaporative  heat  transfer 
in  an  MHP  can  be  provided.  The  liquid  flow  in  a  thin  evaporating  liquid  film  toward 
the  point  x  =  0,  as  shown  in  Fig.  2.3,  takes  place  under  the  influence  of  surface  ten¬ 
sion  and  the  solid-liquid- vapor  interfacial  force  field  (disjoining  pressure).  The  heat 
flux  due  to  conduction  through  the  liquid  in  the  microlayer  region  of  the  meniscus  is 
approximated  as  (Stephan  and  Busse,  1992): 


9  = 


6_ 

Ike 


+ 


Tv(2irRgTvf 5  (2  -  a) 


l-i 


hi 


%Pv 


2a 


(Tw  —  Tfj^y) 


(2.30) 
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Figure  2.3:  Cross  section  of  the  end  of  the  extended  meniscus  region  near  a  heated 


where  the  interfacial  thermal  resistance  is  included.  The  disjoining  pressure  is  (Der- 
jaguin,  1955): 


ft  =  -f  (2.31) 

It  is  more  convenient,  however,  to  describe  the  disjoining  pressure  for  extremely  thin 
films  by  the  following  relation  (Solov’yev  and  Kovalev,  1984;  Holm  and  Goplen,  1979): 


Pd  =  ptRgTs,v  In 


(2.32) 


where  a  and  b  are  constants.  Equation  2.32  was  initially  obtained  by  Potash  and 
Wayner  (1972).  For  pure  water  on  quartz  glass,  a  =  1.534  and  b  =  0.0243  (Holm 
and  Goplen,  1979).  In  the  adiabatic  section,  the  thickness  of  the  liquid  film,  which  is 
extremely  small,  can  be  obtained  from  the  relation  Pd  =  —  cr/ Ftm  by  combining  with 


Eq.  2.32. 

It  is  essential  to  describe  the  difference  between  the  temperature  Ts,v  at  the 
vapor  side  of  the  liquid- vapor  interface  and  the  saturation  temperature  Tv  due  to  the 
disjoining  pressure  of  the  film.  The  following  expression  for  the  pressure  drop  at  the 
surface  of  the  film  in  the  interline  region  is  used  by  Solov’yev  and  Kovalev  (1984), 
which  is  consistent  to  that  by  Carey  (1992). 


In 


Ps,v 


\-Psat(Ts,v ) 


Pd 

P(.Pg  Ts,v 


(2.33) 
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Solov’yev  and  Kovalev  (1984)  obtained  closed-form  solutions  for  q  for  three  regions  of 
an  evaporating  film,  as  shown  in  Fig.  2.3  (see  also  Carey,  1992),  which  are  used  in  the 
present  analysis.  Region  1  represents  the  equilibrium  (nonevaporating)  thin  film,  and 
Region  2  (interline)  is  dominated  by  the  disjoining  pressure.  In  Region  3a  (intrinsic 
meniscus),  the  curvature  of  the  film  varies  under  the  influence  of  surface  tension.  In 
Region  3b  (meniscus),  the  curvature  is  constant.  The  solutions  by  Solov’yev  and 
Kovalev  were  obtained  using  the  boundary-layer  approximation  for  a  liquid  flow,  and 
the  relation  between  the  thickness  and  the  curvature  of  the  film.  The  heat  fluxes  per 
unit  groove  length  can  then  be  estimated  as: 


A  Q2 


1.6pvhfg(Tw  -  Ty) 


hfgpeR2gb(S*  —  <j>p)  j 

.  viHTw  -  Tv)  . 


0.5 


(2.34) 


d0  is  the  equilibrium  (nonevaporating)  thin  film  thickness,  which  can  be  obtained  from 
Eqs.  2.32  and  2.33  at  PgtV  =  Pv  and  T&yV  =  Tw,  where  the  subscript  v  denotes  the  bulk 
conditions  of  the  saturated  vapor.  For  pure  water  on  quartz  glass,  d0  =  1.8  x  10-9  m. 
di  is  defined  from  the  relation  aS\  =  1;  for  water,  Si  =  9.28  x  10~9  m. 


AQ3a 


f±2  (Ty,  -  Tv) 
xi  \  j hpC  T  S (x)  j 


dx 


1  <OPVhjfg(Tyj 


T  \(RA°'Z(  hf 9,762  V’25 

V,\TVJ  \Vlkt(Tw-Tv)J 


(2.35) 
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rw/ ^  (: Tw  -  Tv) 

lx2  l/hpC  +  x2/2Rmki 


A  Q3b  =  r  .■•^-^-4—  dx 

JX'l 


=  i.4(rw  -  rv)(/ipc^i?TO) 


0.5 


7T 


f  hpC82  \ 

vw 


. - arctan 

L2  V  Kt 


0.5i 


(2.36) 


where  hpc  is  the  phase-change  heat  transfer  coefficient  at  the  liquid- vapor  interface 
for  the  intensive  evaporation  process: 


hpc  —  3.2  pvhjg 


(2.37) 


where  it  has  been  assumed  a  —  1. 

In  Eqs.  2.35  and  2.36,  82  is  the  radius  of  the  effective  influence  of  surface 
tension  forces.  For  water,  82  «  10~7  m  (Solov’yev  and  Kovalev,  1984).  The  heat 
transfer  coefficient  of  the  wetted  zone  for  a  given  cross  section  is: 


,  AQ2  +  AQza  +  AQzb  sinq  oq\ 

^e,men  rp  rp  p  .  n\  {Z,.oO ) 

Tw  J.  xi  7tm  cos(7  “f"  vj 

It  should  be  mentioned  that  the  third  term  in  Eq.  2.38  is  dominant.  Then,  for  heat 
input  from  a  hot  fluid,  taking  the  dry  zone  into  account,  it  is  possible  to  estimate  the 
mean  heat  transfer  coefficient  in  the  evaporator  as  follows: 


he 


/  r2i?m  cos(7  +  0) 
U  5  sin  7 


h 


e,men 


-1 


(2.39) 
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The  effective  wall  temperature  in  the  evaporator  in  Eqs.  2.34  -  2.38  is  obtained  in  a 
manner  similar  to  that  in  the  condenser  (Eq.  2.26): 


Tw  —  TeXie  +  (Tv  —  TeX}e)~  f— - f  =-)  (2.40) 

hexye  \’t,ex,e  he' 

When  a  heat  input  in  the  evaporator  is  given  by  Q(z ),  the  dry  zone  of  the 
wall  can  be  overheated  (see  Figs.  2.1  and  2.2(a)).  To  estimate  the  overheating,  the 
thermal  conductance  equation  is  considered: 

,  d2T  1  dQ 

w  wds 2  2N(B/2  +  tw)  dz 

where  T  is  the  mean  cross  sectional  wall  temperature  at  each  point  on  s. 

Eq.  2.41,  it  is  assumed  that  the  heat  flux  on  the  wall  at  each  point  on  z 
The  boundary  conditions  are: 


(2.41) 

In  deriving 
is  constant. 


T\  —  T  ■  — 

U=Ia  ~  w’  ds 


=  0 


s=0 


(2.42) 


Solving  Eq.  2.41,  the  overheat  can  be  expressed  as 


ATe,w  —  T\s-o  —  T\s=Ll  = 


LI 


dQ 


2N(B/2  +  tw)kwtw  dz 
Therefore,  the  maximum  thermal  resistance  of  the  evaporator  is: 


(2.43) 
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T  -T 

-*■  w,max  -*■  v 

Ql 


1 


dQ_ 
Jo  dz 


QaLeN 


_ fi _ +  4=- 

L2(R/2  -T  tw)kwtw  Bhe 


dz  (2.44) 


2.6  Liquid  Flow  in  a  Capillary  Groove 

The  problem  of  liquid  flow  in  triangular  grooves  (Eqs.  2.5  and  2.13)  was  solved  in 
dimensionless  form  by  Ayyaswamy  et  al.  (1974)  and  Khrustalev  (1981).  Ayyaswamy 
et  al.  (1974)  neglected  the  effects  of  shear  stresses  at  the  liquid- vapor  interface.  Since 
the  effects  of  shear  stresses  at  a  free  liquid  surface  may  be  important  for  MHPs,  the 
results  of  Khrustalev  (1981)  are  used,  where  this  two-dimensional  boundary-value 
problem  was  solved  by  the  Ritz  method  in  terms  of  R-functions.  The  solution  gives 
the  values  of  the  dimensionless  volume  flux  through  a  triangular  groove,  which  is 
defined  similar  to  that  given  by  Kamotani  (1976b): 


vMF)  =  JJAt 


wt(x,y)nt 


HW'C 


dPf 


+ 


Pe 


dz  2  Ai  dz 


■r-(t 0(At)  +  pig  sintp 


dx  dy 


(2.45) 


where 


F  = - jo - - -  (2.46) 

2W{iJ  +  JktT^A‘)  +  msi^) 

The  variable  1/V  represents  the  hydraulic  resistance  of  the  groove.  Then,  using  Eqs. 
2.2  and  2.45,  the  pressure  gradient  can  be  defined  as: 
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(2.47) 


At  =  Ri  cos(7  +  S)  (sin(7  +  6)  +  — -)  -  (|  -  7  -  #)  (2.49) 

For  7  =  30°  and  relatively  small  values  of  F,  the  results  given  by  Khrustalev  (1981) 
can  be  approximated  as  follows: 

V(8,F)  =  {ao  +  al0  +  a202)eF  (2.50) 

where  Go  =  0.0241,  Gi  =  0.000905,  02  =  8.75  x  10-6,  and  6  is  in  arc  degrees.  This 
approximation  is  in  agreement  with  the  results  of  Ayyaswamy  et  al.  (1974)  for  F  =  0. 

2.7  Numerical  Treatment 

The  four  main  variables  dependent  on  z  to  be  found  are  as  follows: 

Y  =  R~l:  curvature  of  the  free  surface  of  the  meniscus, 
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m  =  wipiAi :  mass  flow  rate  of  liquid,  in  one  corner, 
po  liquid  pressure, 
pv :  vapor  pressure. 

The  following  system  of  four  ordinary  differential  equations  were  solved  numerically 
using  the  fourth-order  Runge-Kutta  procedure: 


dY  _  1 fdpv_  _  dpt 
dz  a\dz  dz 


(2.51) 
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or  for  a  given  heat-load  function  Q(z ) 
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The  boundary  conditions  at  z  =  0  are: 


pvg  sin  v? 


(2.55) 


Y  =  Rm o,  m  =  0,  pv  =  pv o,  pe  =  pv o  - 

■^mO 


(2.56) 
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The  following  expressions  were  also  used:  Eqs.  2.3,  2.10,  2.18  -  2.27,  2.34  -  2.40, 
2.44,  2.46,  2.48,  2.49  and  2.50. 

For  the  case  of  a  known  heat  load,  it  is  convenient  to  use  the  approximation 
for  the  vapor  pressure  distribution  given  by  Faghri  (1989),  instead  of  numerically 
solving  Eq.  2.55: 


P  =  Pv,  o  +  Pv 
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Le  +  La  <  z  <  Lt  —  Lb  (2.59) 


where  po,a  and  po,c  are  the  values  of  the  pressure  at  the  inlets  of  the  adiabatic  and 
condenser  zones,  respectively.  The  agreement  of  the  results  of  the  numerical  solution 
of  Eq.  2.55  and  Eqs.  2.57  -  2.59  is  within  several  percent  of  the  total  vapor  pressure 
drop  for  the  given  MHP  parameters. 
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An  iterative  procedure  was  used  to  obtain  the  maximum  heat  transfer  ca¬ 
pacity,  Qma.x,  and  the  overall  thermal  resistance,  rt,  of  an  MHP.  The  overall  thermal 
resistance  is: 


rt-  re  +  rc  +  rv 


(2.60) 


where  rv  is  the  thermal  resistance  of  the  vapor  flow,  which  is  defined  as 


QaVL 


1  [L* 


i  rLt-Lb 

r—r  Tv  dz 

■^C  -^b  J Le+La 


(2.61) 


The  mass  of  circulating  fluid,  Mt  =  Mt  -  Mb,  which  represents  the  right- 
hand  side  of  Eq.  2.3  minus  Mb,  and  the  meniscus  radius  Rm  —  Y_1  are  obtained  as 
functions  of  several  parameters: 


Mt  =  Mt(Rmo,  Qa ,  e,  Tv,  W-.)  (2.62) 

Rm  =  Rm(z ,  Rm0,  Qa,  e,  Tv,  ...)  (2.63) 

Analyzing  these  functions  and  taking  into  account  Eq.  2.3,  some  conclusions  con¬ 
cerning  rt,  Lb,  and  Qmax  can  be  made.  For  the  distribution  of  the  liquid  in  an  MHP, 
it  is  logical  to  assume  that,  at  the  maximum  heat  load  ( Qa  —  Q max),  Rmo  has  the 
smallest  value  possible,  provided  that  Eqs.  2.51  -  2.55  can  still  be  solved  successfully 
and  Eq.  2.3  is  satisfied.  It  should  be  noted  that  the  variation  of  the  maximum  heat 
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transport  due  to  decreasing  Rm0  is  very  small  when  Rm0  <  0.1  Dh,v-  The  exact  value 
of  Rm, min,  which  is  different  for  each  MHP  design,  can  be  obtained  by  comparison  of 
the  predicted  and  experimental  results  for  Qmax  for  MHP  operation  against  gravity. 
To  obtain  the  value  for  Lb  and  r,  an  iterative  procedure  is  used,  varying  the  effective 
condenser  length  and  Rm o- 

The  numerical  results  for  the  maximum  heat  transfer  capacity  Qmax  are 
obtained  with  the  following  procedure: 

1.  For  a  definite  small  Rm o  and  Q'a,  the  system  of  Eqs.  2.51  -  2.55  is  solved, 
and  the  obtained  values  of  Me  and  Rm(Lt )  are  analyzed.  If  Me  <C  Mt  and 
Rm(Lt )  <  Rm, max,  liquid  may  be  blocking  part  of  the  condenser.  Therefore, 
the  length  of  the  liquid  blocking  zone  Lb  is  estimated  and  the  new  effective 
condenser  length  L'c  =  Lc  —  Lb  is  set. 

2.  The  value  of  the  total  heat  load  is  increased,  Q"a  —  Q'a  +  A Qa,  and  step 
1  is  repeated  until  Rm(Lt)  ~  i?m; max,  and  the  variation  of  the  meniscus 
contact  angle  6  occurs  near  the  beginning  of  the  blocking  zone.  This  total 
heat  load  is 

3.  The  validity  of  Qmax  is  checked  by  decreasing  R'm0.  If  no  change  occurs, 
go  to  step  4.  Otherwise,  repeat  steps  1  and  2. 

4.  If  there  is  a  lack  of  liquid  in  the  MHP  and  Me  =  Mt  is  obtained  in  the 
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first  two  steps  for  Rm(Lt )  <  -Rm,max,  the  corresponding  maximum  value  of 
Qa  is  considered  to  be  <3max,  provided  that  decreasing  Rm o  has  no  effect 
on  the  results. 

2.8  Results  and  Discussion 

To  verify  the  numerical  results  obtained,  the  experimental  data  provided  by  Wu  and 
Peterson  (1991)  were  used.  Therefore,  the  presented  results  refer  to  an  MHP  with 
the  following  characteristics:  Le  =  0.013  m,  La  =  0.031  m,  Lc  =  0.013  m,  B  =  0.0007 
m,  7  =  30°,  <p  —  0,  Mt  =  0.0032  g,  N  =  4,  tw  =  0.00019  m,  a  =  1.  The  working  fluid 
was  water,  and  two  containers,  made  of  copper  and  silver,  were  examined.  The  axial 
heat  distribution,  Q(z),  was  specified  as: 


Q  a.Z  !  -he  > 

0  <  z  <  Le 

Q  a  5 

L e  ^  Z  ^  Le  d~  La 

(2.64) 

Le  +  La  <  z  <  Lt  —  Lb 

and  Eq.  2.53,  instead  of  Eq.  2.52,  was  used  in  the  numerical  procedure. 

Most  of  the  results  of  the  present  analysis  are  for  heat  loads  close  to  the 
performance  limitations.  The  data  in  Figs.  2.4,  2.5  and  2. 7-2. 9  were  obtained  for 
a  vapor  temperature  at  the  evaporator  end  cap  of  Tv o  =  60° C.  The  hydrodynamic 
performance  characteristics  of  the  MHP  with  no  excess  liquid  is  shown  in  Fig.  2.4. 
Dashed  lines  correspond  to  the  case  of  neglecting  shear  stresses  in  the  liquid  at  the 
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Figure  2.4:  Comparison  of  the  performance  of  the  micro  heat  pipe  with  and  without 
shear  stress  at  the  vapor-liquid  interface  for  the  case  of  no  excess  liquid  in  the  MHP 
(6*o  =  10°)  (a)  liquid  cross  sectional  area,  (b)  radius  of  curvature  of  the  meniscus,  (c) 
pressure  variation. 
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liquid- vapor  interface,  and  solid  lines  refer  to  the  case  of  including  the  effects  of  shear 
stresses  (in  both  cases  the  vapor  distribution  is  the  same).  One  of  the  distinguishing 
features  of  an  MHP  is  the  variable  liquid  cross  sectional  area  along  the  axis.  As  seen 
in  Fig.  2.4(a)  for  60  =  10°,  A(  increased  by  more  than  ten  times  along  the  length 
of  the  MHP,  while  the  radius  of  the  meniscus,  Rm,  changed  by  only  four  times,  as 
shown  in  Fig.  2.4(6).  For  an  inclination  angle  of  <p  =  0,  the  meniscus  curvature 
increased  monotonically  along  the  2  direction.  The  pressure  drop  in  the  liquid  was 
several  times  higher  than  that  of  the  vapor,  which  is  illustrated  in  Fig.  2.4(c).  The 
liquid  pressure  increased  drastically  in  the  evaporator,  where  the  cross  sectional  area 
of  the  liquid  was  small.  Thus,  extending  the  cross  sectional  area  of  the  liquid,  due  to 
the  decrease  of  that  of  the  vapor,  may  increase  the  maximum  heat  transfer  capacity  of 
the  MHP.  The  results  in  Fig.  2.4  are  given  to  show  some  basic  information  concerning 
the  hydrodynamic  characteristics  of  the  MHP  when  the  radius  of  curvature  Rm0  is 
held  constant  at  50  /rm.  It  should  be  noted  that,  in  this  example,  the  mass  of  liquid 
within  the  heat  pipe  is  not  fixed,  but  the  value  of  Rmo  was  the  same  for  all  curves, 
which  gave  the  opportunity  to  show  the  influence  of  the  vapor-liquid  interaction  on 
the  hydraulic  MHP  characteristics  more  clearly. 

The  influence  of  the  liquid  charge  on  Qmax  in  the  case  when  there  is  a  lack 
of  liquid  in  the  MHP  is  illustrated  in  Fig.  2.5  (0O  =  10°).  For  the  considered  MHP, 
this  situation  only  occurs  in  the  absence  of  shear  stress  at  the  free  liquid  surface 
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Figure  2.5:  Dependence  of  the  maximum  heat  transfer  on  (a)  meniscus  contact  angle 
at  the  evaporator  end  cap,  (b)  mass  of  circulating  liquid  in  the  case  of  insufficient 
liquid  fill  (60  =  10°,  Rm0  =  50/zm) 
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and  for  small  90.  The  minimum  contact  angle,  #o,  has  a  significant  influence  on  the 
performance  of  the  MHP,  especially  on  the  maximum  heat  transfer  capacity,  as  shown 
in  Fig.  2.5  for  small  liquid  charges.  According  to  Stepanov  et  al.  (1977),  #0  =  33°  for 
copper-water  and  38°  for  silver-water  pairs,  and  Wu  and  Peterson  (1991)  gave  these 
values  to  be  55°  and  45°,  respectively. 

The  comparison  of  the  numerical  results  and  the  experimental  data  reported 
by  Wu  and  Peterson  (1991)  for  the  maximum  heat  transfer  capacity  are  shown  in  Fig. 
2.6  for  MHPs  with  copper  and  silver  casings.  While  the  data  for  the  minimum  con¬ 
tact  angles  are  contradictory  and  can  be  influenced  by  many  physical  factors,  the 
numerical  results  are  presented  for  9q  =  33°,  38°  and  45°.  The  agreement  between 
the  experimental  results  of  the  onset  of  dryout  and  the  present  numerical  prediction  is 
good.  Neglecting  the  shear  stress  at  the  free  surface  of  the  liquid  due  to  vapor/liquid 
frictional  interaction  (”no  shear  stress”  in  Figs.  2.6  and  2.8)  can  lead  to  an  overes¬ 
timation  of  the  maximum  heat  transfer  capacity.  The  shear  stress  at  the  free  liquid 
surface  influences  the  liquid  distribution  along  the  heat  pipe,  which  can  result  in  an 
increase  of  the  liquid  blocking  of  the  condenser  end  in  comparison  with  the  case  of 
neglecting  this  shear  stress.  The  liquid  distribution  in  an  MHP  is  affected  by  the  heat 
load,  as  shown  in  Fig.  2.7(a)  (#0  =  33°).  The  liquid  distribution  and  the  heat  load 
determine  the  temperature  variation  along  a  heat  pipe  (Fig.  2.7(b)). 

As  seen  in  Fig.  2.8  (  60  —  33°),  the  overall  thermal  resistance  of  the  MHP 
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Figure  2.6:  Maximum  heat  transfer  versus  operating  temperature  (a)  copper-water 
MHP,  (b)  silver- water  MHP 
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Figure  2.7:  Performance  characteristics  of  the  copper-water  MHP  (a)  radius  of  the 
liquid  surface  curvature,  (b)  temperature  distribution 


Figure  2.8:  Influence  of  total  heat  i 
length  (copper-water  MHP) 
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is  summarized  by  three  terms:  re,  rc,  rv,  which  must  be  taken  into  consideration, 
although  the  resistance  of  the  vapor  flow  is  the  greatest  for  the  conditions  studied. 
It  should  be  noted  that,  according  to  experimental  data  by  Wu  and  Peterson  (1991), 
the  thermal  resistance  of  the  MHP  which  was  not  charged  with  working  fluid  (empty 
tube)  was  approximately  150  K/W. 

The  variation  of  the  meniscus  contact  angle  in  the  condenser  is  shown  in 
Fig.  2.9.  In  this  region,  where  progressive  variation  of  the  contact  angle  occurred, 
the  liquid  cross  sectional  height  H  increased  sharply.  Therefore,  the  second  radius  of 
curvature  of  the  liquid  surface,  f?2,  becomes  significant  in  this  region.  In  the  present 
paper,  its  influence  on  the  distribution  of  liquid  along  the  MHP  was  neglected  for 
all  zones.  However,  to  complete  the  physical  model,  it  can  be  assumed  that  the  free 
surface  of  the  liquid  blocking  zone  has  a  curvature  with  a  main  radius  of  Rb,  which 
should  be  estimated  as: 


Rb  "  + 


(2.65) 


Rm  S  z=Lt—  Lb 

From  the  results  of  the  MHP  mathematical  model,  a  new  flat  MHP  configu¬ 
ration  with  an  increased  liquid  cross  sectional  area  for  electronic  cooling  applications 
is  proposed,  as  shown  in  Fig.  2.10.  Due  to  surface  tension  forces, 
liquid  is  contained  in  the  narrow  passages  separating  the  vapor  spaces,  and  is  also 
distributed  along  the  flat  heat-loaded  walls.  The  advantages  of  this  configuration  are 
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Figure  2.10:  New  configuration  cross  section  of  a  flat  MHP 
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an  increase  in  the  maximum  heat  transfer  capacity  and  a  lower  manufacturing  price 
in  compared  to  several  MHPs  investigated  by  Wu  and  Peterson  (1991). 


2.9  Conclusions 

The  numerical  results  of  the  mathematical  model  describing  the  fluid  flow  and  heat 
transfer  within  a  micro  heat  pipe  are  summarized  as  follows: 

1.  Shear  stresses  in  the  liquid  at  the  liquid- vapor  interface  due  to  frictional 
interaction  significantly  influenced  the  maximum  heat  transfer  capacity. 

2.  Accounting  for  shear  stresses  at  the  liquid  free  surface  increased  the 
length  of  the  liquid  blocking  zone  in  the  condenser. 

3.  The  dynamic  component  of  the  pressure  gradient  in  the  liquid  had  no 
pronounced  effect  on  the  performance  characteristics  of  the  MHP. 

4.  For  nearly  maximum  heat  loads,  the  largest  portion  of  the  liquid  pressure 
drop  occurred  in  the  evaporator  and  the  beginning  of  the  adiabatic  section,  where 
the  liquid  cross  sectional  area  was  several  times  smaller  than  that  in  the  condenser. 

5.  The  dominant  thermal  resistances  within  the  MHP  were  those  of  the 
vapor  flow  and  the  liquid  film  in  the  evaporator  and  condenser. 

6.  The  amount  of  working  fluid  and  the  minimum  wetting  contact  angle 
strongly  influenced  the  performance  characteristics  of  the  MHP. 
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7.  The  variation  of  the  meniscus  contact  angle  in  the  condenser  occurred 


near  the  liquid  blocking  zone. 


73 


Chapter  3 


EVAPORATION  AND 
CONDENSATION  ON 
GROOVED  STRUCTURES 

3.1  Summary 

A  detailed  mathematical  model  is  developed  which  describes  heat  transfer  through 
thin  liquid  films  in  the  evaporator  and  condenser  of  heat  pipes  with  capillary  grooves. 
The  model  accounts  for  the  effects  of  interfacial  thermal  resistance,  disjoining  pressure 
and  surface  roughness  for  a  given  meniscus  contact  angle.  The  free  surface  temper¬ 
ature  of  the  liquid  film  is  determined  using  the  extended  Kelvin  equation  and  the 
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expression  for  interfacial  resistance  given  by  the  kinetic  theory.  The  numerical  results 
obtained  are  compared  to  existing  experimental  data.  The  importance  of  the  surface 
roughness  and  interfacial  thermal  resistance  in  predicting  the  heat  transfer  coefficient 
in  the  grooved  evaporator  is  demonstrated. 

3.2  Introduction 

In  an  axially-grooved  heat  pipe  (AGHP),  condensate  return  is  accomplished  by  the 
flow  of  liquid  in  the  grooves  under  the  influence  of  surface  tension,  and  the  maxi¬ 
mum  heat  flux  in  the  evaporator  is  restricted  by  the  fluid  transport  limit  or  by  the 
boiling  limitation.  The  value  of  the  evaporative  heat  transfer  coefficient  is  related 
to  the  maximum  heat  flux  by  the  superheat  of  the  evaporating  liquid  film,  which  is 
critical  for  the  onset  of  nucleate  boiling.  Therefore,  a  detailed  mathematical  model 
is  developed  which  includes  both  the  heat  transfer  through  the  thin  liquid  films  and 
heat  conduction  in  the  fin  between  grooves  and  in  the  meniscus  region  film.  With  this 
model,  the  local  heat  transfer  coefficients  in  the  evaporator  and  condenser  sections  of 
low-temperature  AGHPs  are  determined.  The  present  analysis,  which  incorporates 
several  one-dimensional  boundary-value  problems,  has  the  following  comparatively 
new  features: 
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•  The  heat  transfer  through  the  liquid  films  in  both  the  evaporator  and  condenser 
is  described  with  respect  to  the  disjoining  pressure,  interfacial  thermal  resistance 
and  surface  roughness  or  curvature  (Figs.  3.1  and  3.2). 

•  The  free  surface  temperature  of  the  liquid  film  is  determined  using  the  extended 
Kelvin  equation  and  the  expression  for  interfacial  resistance  given  by  the  kinetic 
theory. 

•  Heat  conduction  in  the  fin  between  grooves  and  the  liquid  meniscus  is  consid¬ 
ered  for  both  the  evaporator  and  condenser  sections  using  a  one-dimensional 
approximation. 

Because  of  these  features  the  proposed  model  is  a  significant  contribution 
over  the  previous  investigators’  attempts  (Kamotani,  1976b,  1978;  Vasiliev  et  ah, 
1981;  and  Stephan  and  Busse,  1992).  Also,  the  model  of  ultra-thin  film  evaporation 
is  employed  which  is  consistent  with  that  given  by  Carey  (1992)  and  Wayner  et  al. 
(1976).  The  present  model  is  developed  for  rectangular,  triangular  and  trapezoidal 
grooves  in  a  circular  tube,  but  flat  evaporators  and  condensers  are  also  described 
by  the  presented  equations.  Heat  transfer  processes  in  the  heat  pipe  container  and 
working  fluid  were  considered  to  be  one-dimensional  in  the  radial  direction,  such 
that  axial  heat  conduction  was  neglected.  The  emphasis  has  been  placed  on  the 
formation  of  the  thin  liquid  films  affected  by  the  operational  conditions.  During  the 
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0 


equilibrium  film 


Figure  3.2:  Thin  evaporating  film  on  a  fragment  of  the  rough  solid  surface 


condensation  process,  liquid  in  the  subcooled  thin  film  flows  towards  the  meniscus 
region  along  the  s-coordinate,  as  shown  in  Fig.  3.1(a).  During  evaporation,  liquid  in 
the  superheated  thin  film  flows  from  the  meniscus  region  in  the  opposite  direction, 
as  presented  in  Figs.  3.1(b)  and  3.2.  The  numerical  results  were  obtained  using 
an  iterative  mathematical  procedure  which  involved  the  following  boundary-value 
problems: 

1.  Formation  of  and  heat  transfer  in  thin  films. 

2.  Heat  transfer  in  the  evaporating  film  on  a  rough  surface. 

3.  Heat  transfer  in  the  condensate  film  on  the  fin  top  surface. 

4.  Heat  conduction  in  a  metallic  fin  and  liquid  meniscus. 

These  problems  are  described  in  detail  in  the  following  sections. 

3.3  Formation  of  and  Heat  Transfer  in  Thin  Liq¬ 
uid  Films 

The  thermal  resistance  of  a  low-temperature  AGHP  depends  mostly  on  the  thickness 
of  the  thin  films  in  the  condenser  and  evaporator  sections.  Since  the  heat  transfer  and 
fluid  dynamics  processes  in  a  thin  film  are  similar  in  both  sections,  it  is  possible  to 
describe  the  formation  of  the  films  by  the  same  equations,  but  taking  into  account  the 


79 


different  directions  of  the  temperature  potential.  In  this  section  a  thin  evaporating 
film  on  a  heat-loaded  surface  with  curvature  I\w  is  considered,  as  shown  in  Fig.  3.2(b). 
The  local  heat  flux  through  the  film  due  to  heat  conduction  is 

Tw-Ts  t  ^ 

q  =  ki - 6 -  l3-1) 

where  the  local  thickness  of  the  liquid  layer  8  and  the  temperature  of  the  free  liquid 
film  surface  Ts  are  functions  of  the  s-coordinate.  For  small  Reynolds  numbers,  an 
assumption  of  a  fully  developed  laminar  liquid  flow  velocity  profile  is  valid: 


dpt  _  _  dK_  dp£  _  dcr  dTs  d_  2  2  \  ( }_  }_ 

ds  ~  ~a  ds  +  ds  ~  dTs  ds  +  Ts^PvVv's)  \ 


(3.3) 


K  is  the  local  interface  curvature,  pd  is  the  disjoining  pressure  (Derjaguin,  1955)  and 
the  last  term  is  the  kinetic  reaction  of  the  evaporating  fluid  pressure.  The  impact  of 
the  last  two  terms  on  the  results  was  found  to  be  negligible  in  the  present  analysis, 
therefore  they  are  omitted  in  the  following  equations. 
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The  continuity  equation  for  the  evaporating  liquid  layer  is 


d[s  q 

—  ut  dr]  =  - - 

as  Jo  ajgpi 


(3.4) 


Substituting  equations  (3.1)— (3.3)  into  equation  (3.4)  gives  the  following  relation  for 
the  thickness  of  the  evaporating  film,  £(s): 


3 


kj(Tw  —  T$) 

hfgprf 


(3.5) 


The  film  surface  curvature  K  is  expressed  in  terms  of  the  solid  surface  curvature  Ku 
and  film  thickness  as 


I<  =  I\w  + 


ds 2  L 


1  + 


^Vl 

is) 


-3/2 


(3.6) 


Following  Potash  and  Wayner  (1972),  a  power-law  dependence  of  pd  on  6  is  given  for 
nonpolar  liquids. 


PJ  =  -A'S-b  (3.7) 

For  water,  however,  the  logarithmic  dependence  is  preferable  (Holm  and  Goplen, 
1979). 

It  is  assumed  that  the  absolute  value  of  the  vapor  core  pressure  at  any  z- 
location  along  the  groove  is  related  to  vapor  temperature  by  the  saturation  conditions 
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Pv  —  Psa.t{^v) 


(3.8) 


and  therefore  can  be  defined  for  a  given  Tv  using  the  saturation  tables. 

The  temperature  of  the  interface  Ts  is  affected  by  the  disjoining  and  capillary 
pressures,  and  also  depends  on  the  value  of  the  interfacial  resistance,  which  is  defined 
for  the  case  of  a  comparatively  small  heat  flux  by  the  following  relation. 


9  =  - 


2<n  \  Ilfg  f  pv  (Psat)i5 


(3.9) 


3- a)  j2xR,  (v/I;  VTi 

pv  and  (psa.t)s  are  the  saturation  pressures  corresponding  to  Tv  in  the  bulk  vapor  and 

at  the  thin  liquid  film  interface,  respectively. 

While  equation  (3.9)  is  used  in  the  present  analysis,  it  seems  useful  to  men¬ 
tion  that  for  the  case  of  extremely  high  heat  fluxes  during  intensive  evaporation  in 
thin  films,  Solov’ev  and  Kovalev  (1984)  have  approximated  the  interfacial  heat  flux 
by  the  following  expression: 


q  =  3.2 yj RgTv[(pS!it)s  ~  Pv]  (3.10) 

Equation  (3.10)  was  derived  with  the  assumption  that  the  accommodation  coefficient 
a  =  1  from  the  expressions  given  by  Labuntsov  and  Krukov  (1977). 

The  relation  between  the  vapor  pressure  over  the  thin  evaporating  film, 
(Psat ) 5 1  affected  by  the  disjoining  pressure,  and  the  saturation  pressure  corresponding 
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to  Ts,  pSat{Ts),  is  given  by  the  extended  Kelvin  equation  (Carey,  1992,  p.  323): 

/  \  frr  \  ["(Psat)(5  —  Pszit(T$)  3 ~  Pd  ~  CrK~\  f  N 

(Psat)«  =  Psat(T5)  exp  ^ - —  — - J  l3-11) 

Equation  (3.11)  reflects  the  fact  that  under  the  influence  of  the  disjoining  and  capillary 
pressures,  the  liquid  free  surface  saturation  pressure  (pSat),$  is  different  from  the  normal 
saturation  pressure  pSzt(Ts)  and  varies  along  the  thin  film  (or  s-coordinate),  while  pv 
and  Tv  are  the  same  for  any  value  of  s  at  a  given  ^-location.  This  is  also  due  to  the 
fact  that  Ts  changes  along  s. 

For  a  thin  evaporating  film,  the  difference  between  ( psat)s  given  by  equation 
(3.11)  and  that  for  a  given  Ts  using  the  saturation  curve  table  is  larger.  This  difference 
is  the  reason  for  the  existence  of  the  thin  nonevaporating  superheated  film,  which  is 
in  the  equilibrium  state  in  spite  of  the  fact  that  Ts  >  Tv. 

Under  steady-state  conditions  the  right-hand  sides  of  equations  (3.1)  and 
(3.9)  can  be  equated. 


T$  =  Tw  + 


$  (  \  hfg  1"  pv  (Ps at) 


h  V2  -a)  J2irRg  \-y/%  VTS 


(3.12) 


Equations  (3.11)  and  (3.12)  determine  the  interfacial  temperature,  Ts,  and  pressure, 
(psat)i5-  Tw  has  to  be  provided  as  an  input  to  the  solution  procedure,  resulting  from 
the  solution  of  the  heat  conduction  problem  in  the  fin  between  the  grooves.  The 


four  boundary  conditions  for  equations  (3.5)  and  (3.6)  must  be  developed  taking  the 
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physical  situation  into  account,  as  shown  in  the  following  sections. 

As  the  liquid  film  thins,  the  disjoining  pressure,  pd,  and  the  interfacial  tem¬ 
perature,  Ts,  increase.  Under  specific  conditions,  a  nonevaporating  film  thickness  is 
present  which  gives  the  equality  of  the  liquid-vapor  interface  and  the  solid  surface 
temperatures,  T$  =  Tw.  This  is  the  thickness  of  the  equilibrium  nonevaporating  film 
S0,  which  can  be  determined  from  equations  (3.11)  and  (3.12).  For  a  nonevaporating 
equilibrium  film  ( q  =  0),  it  follows  from  equation  (3.12)  that 

(Psat)i  =  PvJt^  (3.13) 

Substitution  of  equations  (3.7)  and  (3.13)  into  equation  (3.11)  gives 


So  — 


Psat{TVJ')  PiRgTw\n 


(3.14) 


For  water,  the  following  equation  for  the  disjoining  pressure  was  used  (Holm  and 
Goplen,  1979) 


pd  =  piRgTs  In 


3.3 


(3.15) 


where  a  =  1.5336  and  b  =  0.0243.  The  thickness  of  the  equilibrium  film  is  given  for 
water  by 
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dn  =  3.3-! 


1  f^sat  (Tw)-pvJTJTv+aK 

-exp  - DT - +ln 

Cl  Pt  -ILg  -L  yj 


Pv 


Psat(7u,) 


1/b 


(3.16) 


3.4  Heat  Transfer  in  the  Thin-Film  Region  of  the 
Evaporator 

This  problem  has  been  treated  numerically  and  experimentally  by  different  authors, 
whose  results  are  mentioned  here  to  understand  the  basis  of  the  present  model. 
Kamotani  (1978),  Holm  and  Goplen  (1979),  and  Stephan  and  Busse  (1992)  mod¬ 
eled  an  evaporating  extended  meniscus  in  a  capillary  groove  (Fig.  3.1).  In  all  of 
the  above  papers,  it  is  emphasized  that  most  of  the  heat  is  transferred  through  the 
region  where  the  thickness  of  the  liquid  layer  is  extremely  small.  The  significance  of 
the  temperature  difference  between  the  saturated  vapor  core  and  the  interface  has 
been  stressed  by  Solov’yev  and  Kovalev  (1984)  and  Stephan  and  Busse  (1992).  In 
the  mathematical  models  of  the  above  authors,  the  solid  surface  was  assumed  to  be 
smooth.  Kamotani  (1978)  noted  that  this  assumption  and,  generally  including  of  the 
heat  transfer  in  the  microfilm  region  in  the  model  could  lead  to  an  overestimation  of 
the  total  heat  transfer  coefficient.  The  same  concern  has  been  expressed  by  Vasiliev 
et  al.  (1981). 

In  light  of  these  findings,  in  the  present  analysis  the  difference  between  the 
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saturated  vapor  temperature  and  that  of  the  free  liquid  surface  was  considered,  and 
the  existence  of  the  surface  roughness  and  its  influence  on  evaporative  heat  transfer 
was  taken  into  consideration.  In  general,  manufacturing  processes  always  leave  some 
degree  of  roughness  on  the  metallic  surface.  Alloys  of  copper,  brass,  steel  and  alu¬ 
minum  invariably  have  some  distinct  grain  structure,  resulting  from  processing  the 
materials.  In  addition,  corrosion  and  deposition  of  some  substances  on  the  surface 
can  influence  its  microrelief.  This  means  the  solid  surface  is  totally  covered  with  mi¬ 
croroughnesses,  where  the  characteristic  size  may  vary  from,  for  example,  Rr  —  10~8 
to  10"6  m.  Apparently,  the  thin  liquid  film  formation  can  be  affected  by  some  of  these 
microroughnesses.  It  can  be  assumed  that  at  least  some  part  of  a  single  roughness 
fragment,  on  which  the  thin  film  formation  takes  place,  has  a  circular  cross  section 
and  is  extended  in  the  2-direction  due  to  manufacturing  the  axial  grooves  (Fig.  3.2). 

In  the  present  analysis  the  free  liquid  surface  is  divided  into  four  regions 
(Fig.  3.2).  The  first  region  is  the  equilibrium  nonevaporating  film.  The  second 
(microfilm)  region  ranges  in  the  interval  8q  <  8  <  Si,  where  the  increase  of  the  liquid 
film  thickness  up  to  the  value  <fi  is  described  by  equations  (3.5)  and  (3.6).  In  this 
region,  the  generalized  capillary  pressure  pcap  =  aK  —  (here  pcap  was  defined  so 
that  its  value  is  positive)  is  changing  drastically  along  the  s-coordinate  from  the  initial 
value  up  to  an  almost  constant  value  at  point  si,  where  the  film  thickness,  8i,  is  large 
enough  to  neglect  the  capillary  pressure  gradient.  It  is  useful  to  mention  that  some 
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investigators  have  denoted  this  microfilm  region  as  the  “interline  region.”  The  third 
(transition)  region,  where  the  liquid- vapor  interface  curvature  is  constant,  is  bounded 
by  <$i  <  6  <  Rr  +  So,  and  the  local  film  thickness  is  determined  by  the  geometry  of  the 
solid  surface  relief  and  the  value  of  the  meniscus  radius  Rm.  In  the  fourth  (meniscus) 
region,  where  by  definition  S  >  Rr  +  60,  the  local  film  thickness  can  be  considered 
independent  of  the  solid  surface  microrelief.  In  the  third  and  fourth  regions,  the  heat 
transfer  is  determined  by  heat  conduction  in  the  meniscus  liquid  film  and  the  metallic 
fin  between  the  grooves.  However,  in  the  second  region,  the  temperature  gradient  in 
the  solid  body  can  be  neglected  in  comparison  to  that  in  liquid  due  to  the  extremely 
small  size  of  this  region. 

The  total  heat  flow  rate  per  unit  groove  length  in  the  microfilm  region  is 
defined  as 


Equations  (3.5)— (3.8) ,  (3.11)  and  (3.12)  must  be  solved  for  four  variables:  6,  S',  pcap 
and  Q{njc(s)  in  the  interval  from  s  =  0  to  the  point  s  =  Si,  where  pcap  can  be 
considered  to  be  constant.  Now,  instead  of  the  two  second-order  equations  (3.5) 
and  (3.6),  the  following  four  first-order  equations  should  be  considered  with  their 
respective  boundary  conditions: 
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(3.18) 


dS_ 

ds 


=  S' 


cW 

ds 


(3.19) 


d-Pc ap  _  p/  /  \ 

ds  _  hJg6^™c[  ’ 


(3.20) 


dQ'nac  _  ^ 

ds  S/k( 


(3.21) 


<5U=0  =  d0 


(3.22) 


<a=0  =  o 


(3.23) 


Pcap  |  s=C 


Rr  +  £0 


+  axb 


(3.24) 


<4i=U  =  0  (3.25) 

The  value  of  <50  is  found  from  equation  (3.14),  where  I\  =  — 1  /Rr. 

Though  the  initial-value  problem,  equations  (3.18)— (3.25)  is  completely  de¬ 
termined,  its  solution  must  satisfy  one  more  condition: 
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Since  the  only  parameter  which  is  not  fixed  in  this  problem  is  connected  with  the 
surface  roughness  characteristics,  the  boundary  condition  (3.26)  can  be  satisfied  by 
the  choice  of  Rr.  Physically,  it  means  that  the  beginning  of  the  evaporating  film 
is  shifted  along  the  rough  surface  depending  on  the  situation  so  as  to  satisfy  the 
conservation  laws.  However,  in  a  smooth  surface  model  (Rr  — >  oo)  the  solution  will 
probably  not  satisfy  equation  (3.26).  As  a  result  of  this  problem,  the  values  of  and 
can  be  determined  and  the  transition  region  can  be  considered,  provided 
that  61  <  RT,  where  the  free  liquid  surface  curvature  is  constant  and  its  radius  Rm  is 
many  times  larger  than  Rr.  Based  on  the  geometry  shown  in  Fig.  3.2,  the  following 
approximation  for  the  liquid  film  thickness  in  the  transition  region  x j  <  x  <  xtr  is 
given: 


8  =  d0  +  Rr  -  \]Rr2  -  x 2  -  Rm  +  [Rm2  +  x2  +  2 Rmx  sin  0/)  /  (3.27) 

Equation  (3.27)  is  valid  for  the  rough  surface  model  ( Of  can  be  set  equal  to  zero  for 
very  small  Rr)  and  also  the  smooth  surface  model  in  the  meniscus  region  (Rr  — *•  oo 
and  Of  is  given  as  a  result  of  the  microfilm  problem  solution). 

The  heat  flow  rate  per  unit  groove  length  in  the  transition  region  is 


where  x j  and  xtr  are  obtained  from  equation  (3.27)  provided  ^  =  <5i  and  6  —  RT  +  <50, 
respectively. 

Now  the  connecting  point  between  the  transition  and  meniscus  regions  must 
be  considered.  At  this  point,  the  film  thickness,  the  free  surface  curvature,  and  the 
liquid  surface  slope  angle  must  coincide  from  both  sides.  In  the  rough  surface  model, 
the  last  condition  is  always  satisfied  because  the  length  of  the  microfilm  region  is 
smaller  than  Rr,  and  the  rough  fragment  with  the  film  can  be  “turned”  around  its 
center  in  the  needed  direction  (see  Fig.  3.2).  In  other  words,  because  of  the  circular 
geometry  of  the  rough  fragment  and  the  constant  temperature  of  the  solid  surface  in 
the  microfilm  region,  the  slope  of  the  film  free  surface  is  not  fixed  in  the  mathematical 
model.  On  the  contrary,  in  the  smooth  surface  model  the  numerical  results  give  Of 
which  is  generally  not  equal  to  #men  determined  by  the  fluid  flow  along  the  groove. 
Stephan  (1992)  seems  to  have  answered  this  contradiction  using  a  rounded  fin  corner, 
however,  this  explanation  is  not  completely  satisfactory.  Note  that  in  the  situation 
when  Of  ^  0men,  the  smooth  surface  model  can  be  used  along  with  the  rounded  fin 
corner,  where  the  radius  is  Rfa.  In  this  case  equation  (3.27)  can  also  be  used  provided 
Rr  is  changed  to  Rf^. 


It  is  useful  to  mention  here  that  the  values  of  Rm  and  dmen  are  connected 


by  the  geometric  relation  #men  =  arccos  (W/ Rm)  —  7  and  should  be  given  as  a  result 
of  the  solution  of  the  problem  for  the  fluid  transport  along  the  groove.  The  fin  top 
temperature  Tw  should  be  defined  from  the  consideration  of  the  heat  conduction 
problem  in  the  fin  between  grooves  and  in  the  meniscus  liquid  film  discussed  below. 

Simplified  model  of  heat  transfer  in  the  evaporating  thin  film.  The  free 
liquid  surface  curvature  K  in  the  microfilm  region  varies  from  the  initial  value  to 
that  in  the  meniscus  region.  Its  variation  is  described  by  equations  (3.18)— (3.26)  with 
respect  to  the  pcap  and  pd  definitions.  In  spite  of  a  sharp  maximum  which  the  K 
function  has  in  the  microfilm  region,  its  variation  only  slightly  affects  the  total  heat 
transfer  coefficient.  To  check  this  hypothesis  numerically,  a  simplified  version  of  the 
heat  transfer  model  of  the  microfilm  region  was  developed,  where  it  was  assumed 
that  the  microfilm  free  surface  curvature  is  equal  to  that  in  the  meniscus  region. 
Therefore,  instead  of  solving  equations  (3.18)— (3.26),  the  microfilm  thickness  in  this 
region  (and  also  in  the  transition  region)  can  be  given  by  equation  (3.27)  for  the 
interval  0  <  x  <  xtT.  In  this  case,  the  heat  flow  rate  per  unit  groove  length  in  both 
the  microfilm  and  transition  regions  is 

Pnijjr*  <3-29> 
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3.5  Heat  Transfer  in  the  Thin  Film  Region  of  the 


Condenser 

Heat  transfer  during  condensation  on  a  grooved  surface  has  been  considered  by 
Kamotani  (1976b)  and  Babenko  et  al.  (1981)  for  the  case  of  rectangular  and  trape¬ 
zoidal  fins  with  rounded  corners,  respectively,  with  the  assumption  that  Ts  =  Tv. 
Analyzing  their  results,  the  following  conclusions  are  made,  which  lead  to  the  simpli¬ 
fication  of  equations  (3.5)  and  (3.6): 

1.  The  surface  of  the  liquid  film  is  smooth  and  the  film  thickness  variation 
along  s-coordinate  is  weak  (see  Fig.  3.1): 


'dsy 

Ts)  «x 


2.  The  disjoining  pressure  gradient  along  the  film  flow  can  be  neglected  in 
comparison  to  that  of  the  capillary  pressure  due  to  the  surface  tension 
force  because  of  the  large  film  thickness. 


Taking  the  above  points  into  consideration,  and  substituting  equation  (3.6)  into  equa¬ 
tion  (3.5)  gives  the  following  differential  equation  for  the  film  thickness  at  the  top  of 
the  fin  between  grooves: 
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4l<  +  1r)]  = 


3  nikt 
a  hjgpt 


(Ts  -  Tw) 


(3.30) 


The  boundary  conditions  for  equation  (3.30)  at  s  —  0  are 


d3S  n 
0  ’  is3  ~ 


(3.31) 


These  conditions  imply  that  the  thickness  and  curvature  of  the  film  are  symmetric 
around  5  =  0.  For  small  «  (see  Fig.  3.1(a))  at  s  =  C2,  the  curvature  of  the  film  and 
its  surface  slope  angle  are  determined  by  the  radius  of  the  meniscus  in  the  groove: 


ds 2  Rm 


(3.32) 


d8  (  .  W  .  L\ 

—  =  tan  k  —  arcsm  — - h  arcsm  — 

ds  \  ALjji  diy 


(3.33) 


where  L2  is  the  length  of  the  film,  which  is  equal  to  Li  in  the  case  of  a  flat  fin  top 
geometry  (Fig.  3.1).  The  boundary  value  problem,  equations  (3.30)— (3.33),  is  solved 
approximately  by  introducing  the  following  polynomial  function  for  the  film  thickness: 


S(s)  =  Co  +  C\(s  —  L2)  +  C2(s  —  L2)1  +  Cz(s  —  L2)3  +  C*4(s  —  L2) 


(3.34) 


From  the  boundary  conditions  (3.31)— (3.33)  the  values  of  the  coefficients  are 


r  ,  (  •  w  ,  ■  Li\  r  1 

Ci  =  tan  I  k  —  arcsm  — - h  arcsm  —  I  ,  C2  =  „ 

\  •K'V ' 
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„  2 C2L2  -Ci  „  Cz 

°3  2  L\  ’  4T2 

At  the  point  s  —  L2,  where  the  thickness  of  the  film  is  usually  at  a  minimum,  equation 
(3.30)  must  be  satisfied  exactly,  and  the  total  mass  flow  rate  of  the  condensate  due 
to  the  surface  tension  force  must  be  equal  to  the  total  amount  of  fluid  condensed  in 
the  region  0  <  5  <  L2.  Thus,  integrating  equation  (3.30)  we  have 


a  hfgPe 

3 pi 


S3 


dK„ 

dsz  ds 


l 


Ts  -  Tu 


ds 


(3.35) 


-I  (s=Li) 

Substituting  equation  (3.34)  into  equation  (3.35)  and  solving  numerically  for  Co,  the 


heat  flow  rate  per  unit  groove  length  through  the  thin  film  region  is 


Qs  =  / 

Jo 


h(Ts  —  Tu,) 


Co  +  (7i(s  —  T2)  +  Ci^s  —  L2y  +  (^3(5  —  L2y  +  C^s  —  L2y 


ds  (3.36) 


The  fin  top  temperature  Tw  is  given  from  the  results  of  the  heat  conduction  problem 
in  the  fin  and  the  meniscus  region.  Equation  (3.35)  must  also  be  solved  within  the 
following  iterative  procedure  because  of  the  influence  of  the  film  surface  curvature 
and  the  disjoining  pressure  on  Tg.  In  the  first  iteration,  Ts  is  defined  from  equations 
(3.11)  and  (3.12)  assuming  that  K  =  I\w  and  pd  —  0.  In  the  second  and  following 
steps 


K  —  Kw  + 


dPT 

ds 2 
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where  the  last  term  is  calculated  using  the  solution  of  the  previous  step  for  6(s)  and 
pd  =  pd(S)  (here  the  bar  denotes  an  average  value).  While  the  influence  of  K  on  the 
presented  results  was  negligible  in  comparison  with  the  effect  of  the  meniscus  radius 
variation,  it  can  be  important  for  extremely  thin  films  of  condensate  with  large  free 
surface  curvatures,  in  which  case  the  problem  should  be  treated  numerically  in  the 
frames  of  a  more  complicated  analysis. 

Now,  the  consideration  of  the  meniscus  region  gives  the  opportunity  to  ob¬ 
tain  the  heat  transfer  coefficients. 


3.6  Heat  Conduction  in  the  Metallic  Fin  and  Menis¬ 
cus  Region  Film 

For  low-temperature  heat  pipes,  the  thermal  conductivity  of  the  metallic  casing  is 
several  hundred  times  higher  than  that  of  the  liquid  working  fluid.  Nearly  all  of  the 
heat  is  transferred  from  the  metallic  fin  between  grooves  to  the  saturated  vapor  or  vice 
versa  through  a  thin  liquid  film  in  the  vicinity  of  the  fin  top.  The  temperature  drop  in 
the  metallic  fin  is  many  times  smaller  than  in  the  liquid  film  (Schneider  et  al.,  1976; 
Stephan  and  Busse,  1992).  Therefore,  in  the  present  analysis  the  temperature  gradient 
in  the  metallic  fin  in  the  direction  transverse  to  the  ^-coordinate  is  neglected  (Fig. 
3.1).  The  heat  conduction  in  the  metallic  fin  and  meniscus  liquid  film  is  described 
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by  the  following  equation,  which  was  obtained  as  a  result  of  an  energy  balance  over 
a  differential  element  (Vasiliev  et  ah,  1981): 


<PT  dT  tan(7  +  x)  .  (T  T)  h  _  n 

dx2  dx  Lsnix)  1  ;^d(x)Xfin(x)  1  ; 

The  fin  thickness  variation  is  due  to  its  wall  inclination  angle  and  the  circular  tube 
geometry 


Tfin(x)  =  Lx  +  x  tan(7  +  x) 


(3.38) 


and  the  liquid  film  thickness  is 


6  —  62  —  Rm  + 


Rm  2  + 


+ 


2  RmX 


cos2  (7  +  x)  cos(7  +  x) 


sin  Or, 


1/2 


(3.39) 


where  <5  is  measured  perpendicularly  from  the  liquid-vapor  interface.  Equations 
(3.37)— (3.39)  are  valid  for  the  evaporator  and  condenser  sections.  However,  the 
boundary  conditions  for  equation  (3.37)  in  these  two  sections  are  different,  and  the 
value  of  62  should  be  chosen  as  follows: 

82  =  Rr  +  do  in  the  rough  surface  evaporation  model, 

82  =  8 1  in  the  smooth  surface  evaporation  model, 

^2  =  8  |s=i,2  for  the  condenser  heat  transfer  model. 

The  boundary  conditions  for  equation  (3.37)  in  the  evaporator  are 
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dT 

dx 


r=0 


dT 

dx 


Qmic  F~  Qtr 

kvjT\ 


QeTT  Rq 


(3.40) 


(3.41) 


i*=i,  KNIL-i  +  t,  tan(7  +  x)] 
where  boundary  condition  (3.40)  is  written  with  the  assumption  that  xtT  -C  tg.  For 
the  simplified  model,  equation  (3.37)  was  solved  also  in  the  microfilm  and  transition 
regions,  where  8  was  given  by  equation  (3.27)  and  the  right-hand  side  of  equation 
(3.40)  was  set  equal  to  zero. 

The  boundary  conditions  for  equation  (3.37)  in  the  condenser  are 


dT 

dx 


x=0 


dT 

dx 


Q's 

kwLi 


qcnRo 


(3.42) 


(3.43) 


X=tn 


kwNlLx  +  tg  tan(q  +  x)] 

While  the  values  of  Q !adc,  Q'tr  and  Q's  depend  on  Tw  =  T\x=0,  which  is 
obtained  from  the  solution  of  equations  (3 .37)— (3 .43) ,  this  problem  is  to  be  solved  in 
conjunction  with  those  concerning  heat  transfer  in  the  thin  film  regions. 

The  local  heat  transfer  coefficient  (for  a  given  z)  in  the  evaporator  from  the 
bottom  of  the  groove  surface  to  the  vapor  is 
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(3.44) 


he, hot  — 


Ro 


[T\X=tg  —  Ty ]  Ry  +  tt  +  tg 
The  local  heat  transfer  coefficient  from  the  external  surface  of  the  evaporator  to  the 
vapor  is 


he  = 


Ro 


In 


Ro 


+  =- 


1 


Ro 


i  -i 


ikyj  Ry+ti  +  tg  /le,bot  + 

where  the  thermal  resistance  of  the  circular  tube  wall  is  taken  into  account. 


(3.45) 


For  the  condenser  region,  the  heat  transfer  coefficients  are  defined  in  a  sim¬ 
ilar  manner: 


he,  bot  — 


9c 


R0 


[T|(X  =  ts)  —  Tv]  Rv  +  tt  +  tg 


(3.46) 


he  — 


R0 

iky. 


In 


Ro 


+  = 


1 


Ro 


i-i 


Rv  +tt  +  tg  hC  ,bot  Rv  +  t,  +  tgi 


(3.47) 


3.7  Numerical  Treatment 

Equations  (3.11)  and  (3.12)  were  simultaneously  solved  for  T$  (absolute  error  Aa  = 
0.0001  K)  and  (psat)s  (Aa  =  1  Pa)  for  every  point  on  s  by  means  of  Wegstein’s  iteration 
method  (Lance,  1960).  The  system  of  the  four  first-order  ordinary  differential  equa¬ 
tions  with  four  initial  conditions  and  one  constitutive  condition  describing  the  evap- 
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orating  microfilm  region,  equations  (3. IS)— (3.26),  were  solved  using  the  fourth-order 
Runge-Kutta  procedure  and  the  shooting  method  (on  parameter  Rr).  The  controlled 
relative  error  was  less  than  0.001%  for  each  of  the  variables.  The  results  obtained 
for  comparatively  small  temperature  drops  through  the  thin  film  were  compared  with 
those  from  the  simplified  model.  Since  the  agreement  was  good,  the  simplified  model 
was  used  further  in  the  prediction  of  the  AGHP  characteristics.  Equation  (3.35) 
was  solved  for  Co  by  means  of  Muller’s  iteration  method  (Aa  =  I0~n  m),  and  the 
integration  in  equation  (3.36)  was  made  using  Simpson’s  method.  The  heat  conduc¬ 
tion  problem,  equations  (3.37)— (3.43),  was  also  solved  by  the  standard  Runge-Kutta 
method  (Aa  =  0.0001  K  and  Ar  =  0.001%  for  the  functions  T  and  dT/dx,  respec¬ 
tively)  along  with  equations  (3.35)  and  (3.36)  within  the  iterative  procedure  to  find 
Tw  (Aa  =  0.0001  K). 

3.8  Results  and  Discussion 

To  verify  the  numerical  results  obtained,  the  experimental  data  provided  by  Schlitt 
et  al.  (1974)  were  used.  Therefore,  the  presented  results  mostly  refer  to  the  AGHP 
with  the  following  geometry:  Lt  —  0.914  m,  Lc  =  0.152  m,  0.15  <  Le  <  0.343  m, 
W  =  0.305  mm,  tg  =  1.02  mm,  L\  =  0.215  mm,  Rv  =  4.43  mm,  R0  =  7.95  mm, 
7  =  0,  N  =  27,  ac  =  0,  tt  =  0.  The  working  fluids  were  ammonia  and  ethane,  the 
casing  material  thermal  conductivity  was  assumed  to  be  kw  =  170  W/(m-K),  a  —  1 
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(if  another  value  is  not  indicated  in  the  text),  dispersion  constant  A'  —  10  21  J  and 
B  =  3. 

The  data  in  Figs.  3. 3-3. 5  were  obtained  for  ammonia  with  a  vapor  temper¬ 
ature  in  the  evaporator  of  Tv  =  250  K  and  a  =  1.  The  solid  surface  superheat  is 
AT  =  | Tw  —  Tv  |,  and  the  results  obtained  using  the  simplified  model  for  evaporating 
film  are  denoted  as  SIMPL. 

Figure  3.3(a)  shows  the  variations  of  the  free  liquid  surface  temperature 
along  the  evaporating  film  for  AT  =  0.047  K,  0.070  K  and  0.120  K,  which  are  from 
the  solutions  of  equations  (3.7),  (3. 1 1 )— (3.13) ,  and  (3 . 1 8)— (3 .26)  in  the  microfilm 
region.  These  results  are  compared  to  those  obtained  by  the  simplified  model,  where 
equations  (3.7)  and  (3.11 )— (3.13)  were  solved  along  with  equations  (3.27),  (3.37)  and 
(3.38)  with  the  boundary  conditions 

T\X  =  0  =  Tyj,  =0 

ax  x=o 

in  the  microfilm  and  transition  regions  for  the  same  values  of  the  roughness  charac¬ 
teristic  sizes  (Rt  =  0.33  fim,  1.0  ^m  and  Rr  — ►  oo).  It  should  be  noted  that  the 
temperature  drop  in  the  solid  body  in  these  regions  was  negligible  in  the  results  of 
the  simplified  model  in  comparison  to  AT,  and  the  equilibrium  film  thickness  was 
defined  within  the  assumption  that  its  free  surface  curvature  was  equal  to  1  / Rm.  In 
the  simplified  model  for  the  case  of  a  smooth  surface,  the  value  of  the  contact  angle 
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Figure  3.3:  Characteristics  of  the  evaporating  film  along  the  solid-liquid  interface 
(ammonia,  Tv  =  250  K)  (a)  free  liquid  surface  temperature,  (b)  thickness  of  the  film, 
(c)  generalized  capillary  pressure 
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Figure  3.5:  Local  heat  transfer  coefficient  in  the  evaporator  of  the  ammonia- A1  heat 
pipe  (Tv  =  250  K)  (a)  versus  roughness  size,  (b)  versus  heat  flux  ( Rr  —  1  ^m) 
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in  the  microfilm  region  was  0f  =  7°,  which  was  given  by  the  numerical  solution  of 
equations  (3.18)— (3.26). 

The  corresponding  variations  of  the  film  thickness  8  and  generalized  capillary 
pressure  pcap  are  shown  in  the  Figs.  3.3(b)  and  3.3(c).  The  results  obtained  by  the 
simplified  model  have  been  artificially  shifted  along  the  ^-coordinate  in  these  figures 
(and  also  in  Fig.  3.4(a))  to  make  the  comparison  more  understandable.  Also,  it 
should  be  noted  that  there  is  some  difference  between  the  s-coordinate  and  the  x- 
coordinate  used  in  the  simplified  model.  The  following  relation  has  been  used  in  the 
present  paper:  s  —  RT  arcsin  ( x/Rr ). 

In  Fig.  3.3(a),  the  interval  of  Ts  variation  along  the  evaporating  film  from 
the  value  of  Tw  to  approximately  Tv  was  more  prolonged  in  comparison  to  the  re¬ 
sults  by  Stephan  and  Busse  (1992),  and  the  interfacial  thermal  resistance  was  still 
significant  even  when  the  film  thickness  was  larger  than  0.1  pm.  For  a  smaller  char¬ 
acteristic  size  Rt,  the  film  thickness  increased  more  sharply  along  the  solid  surface 
(Fig.  3.3(b)),  which  is  in  agreement  with  equation  (3.27).  It  should  be  mentioned 
that  for  the  problem,  equations  (3.18)— (3.26)  (unlike  for  the  simplified  model)  Rr  is 
not  a  parameter  but  the  result  of  the  numerical  solution.  The  values  of  the  maximum 
heat  flux  in  the  microfilm  region  were  extremely  high  in  comparison  to  those  in  the 
meniscus  region  (Fig.  3.4).  For  AT  =  0.120  K,  the  generalized  capillary  pressure 
Pcap  decreased  from  the  initial  value  to  the  almost  constant  by  approximately  5000 
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times  (Fig.  3.3(c)).  For  a  larger  AT,  this  sharp  decrease  can  cause  some  difficulties 
in  the  numerical  treatment  while  solving  equations  (3. 18)— (3.26);  that  is  why  the 
simplified  model  is  useful.  The  simplified  model  has  given  the  variation  of  pcap  along 
the  film  which  is  even  more  drastic  because  of  the  surface  tension  term  is  absent 
in  the  capillary  pressure  gradient  (Fig.  3.3(c)).  However,  the  decrease  of  the  total 
heat  flow  rate  in  the  microfilm  region  caused  by  this  assumption  was  comparatively 
small,  which  is  illustrated  by  Fig.  3.4(a).  The  distributions  of  the  heat  flux  in  the 
microfilm,  transition  and  beginning  of  meniscus  regions  for  different  meniscus  contact 
angles  6men  as  predicted  by  the  simplified  model  are  presented  in  Fig.  3.4(b).  The 
total  heat  flow  through  the  meniscus  region  was  significantly  larger  in  comparison 
to  that  through  the  microfilm  region.  This  means  that  while  estimating  the  heat 
transfer  coefficient  for  an  evaporator  element,  shown  in  Fig.  3.1,  the  simplified  model 
should  provide  the  accuracy  needed.  To  verify  this,  the  numerical  results  for  the  local 
heat  transfer  coefficient  he  in  Fig.  3.5(a)  have  been  obtained.  The  simplified  model 
underestimated  he  by  only  5%,  which  enables  its  use  when  it  is  necessary  to  avoid  the 
numerical  difficulties  mentioned  above.  The  local  evaporative  heat  transfer  coefficient 
he  depends  upon  the  meniscus  contact  angle  #men,  especially  for  small  9me n,  and  is 
practically  independent  of  the  heat  flux  on  the  external  wall  surface  of  the  evaporator 
and  also  of  AT,  as  shown  in  Fig.  3.5(b).  The  characteristic  roughness  size  affected 
the  value  of  he,  decreasing  it  up  to  30%  for  a  =  1  in  comparison  to  the  value  obtained 
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for  the  smooth  solid  surface.  For  large  meniscus  contact  angle  the  influence  of  the 
roughness  size  on  the  heat  transfer  coefficient  is  at  the  maximum  when  Rr  is  close  to 
the  length  of  the  microfilm  region.  For  small  values  of  the  accommodation  coefficient 
(for  example  for  a  =  0.05)  the  effect  of  the  surface  roughness  on  the  heat  transfer  is 
insignificant  because  the  heat  flux  in  the  microfilm  region  in  this  case  is  comparatively 
small  (Fig.  3.6). 

The  results  of  the  present  model  were  compared  with  the  experimental  data 
by  Schlitt  et  al.  (1974)  and  Ivanovskii  et  al.  (1984)  for  the  case  of  a  small  heat  load 
applied  to  the  AGHP  (or  evaporator).  For  a  small  heat  load  ( Qa  <C  Qmzx )  the  values 
of  the  meniscus  angle  in  both  the  evaporator  and  condenser  of  the  AGHP  under 
consideration  are  comparatively  large:  9men  >  60°  in  the  evaporator  and  #men  >  80° 
in  the  condenser.  This  is  valid  because  in  the  case  without  a  heat  load  the  grooves 
of  an  AGHP  in  the  horizontal  position  are  completely  filled  with  liquid  (i.e. ,  the 
meniscus  angle  is  close  to  90°).  For  6men  >  60°  the  local  evaporative  heat  transfer 
coefficients  are  practically  independent  on  #men,  as  shown  in  Fig.  3.6.  The  values  of 
the  evaporative  heat  transfer  coefficients  (based  on  the  outer  tube  diameter)  obtained 
experimentally  by  Schlitt  et  al.  (1974)  and  those  reported  by  Ivanovskii  et  al.  (1984) 
were  also  found  to  be  independent  of  heat  load,  which  resulted  in  a  valid  comparison, 
as  given  in  Table  3.1.  The  agreement  of  the  results  for  ammonia,  ethane  and  water 
is  good  for  a  <C  1  since  it  was  mentioned  by  Carey  (1992,  p.  119)  that,  for  some 
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Figure  3.6:  Effect  of  the  meniscus  contact  angle  on  the  local  evaporative  heat  transfer 
coefficients  (AT  =  1  K)  (a)  ammonia-Al  heat  pipe  by  Schlitt  et  al.  (1974),  (Tv  =  250 
K),  (b)  ethane- Al  heat  pipe  by  Schlitt  et  al.  (1974),  ( Tv  =  200  K),  (c)  water-copper 
evaporator  by  Ivanovskii  et  al.  (1984),  (T„  =  300  K) 
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Table  3.1:  Comparison  of  the  predicted  he  with  existing  experimental  data 


Investigators 

Schlitt  et  al.  (1974) 

Schlitt  et  al.  (1974) 

Ivanovskii  et  al.  (1984) 

Working  fluid 

Ammonia 

Ethane 

Water 

Tv  (K) 

250 

200 

300 

Casing  material 

Aluminum 

Aluminum 

Copper 

Experimental  value  of 

he  (W/m2-K) 

3920 

770 

9500 

Present  prediction, 

K  (W/m2-K), 

9  =  60°,  a  =  0.05 

4140 

1180 

8620 
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substances  (ethanol,  methanol,  water,  etc.),  the  accommodation  coefficient  had  been 
found  to  have  very  small  values  (0.02  to  0.04)  in  the  experiments  by  Paul  (1962).  The 
physical  reason  for  low  a  values  in  the  microfilm  region  of  the  evaporator  can  be  the 
concentration  of  the  contaminants  which  usually  exist  in  a  heat  pipe  in  this  region. 
For  the  case  of  a  =  1,  the  prediction  gave  significant  (up  to  100%)  overestimations  of 
he  even  for  a  rough  surface,  as  can  be  seen  from  Fig.  3.6.  The  experimental  data  by 
Ivanovskii  et  al.  (1984)  correspond  to  the  case  of  evaporation  of  water  from  a  copper 
plate  with  rectangular  grooves  for  heat  fluxes  on  the  wall  up  to  20  W/cm2  ( W  =  0.17 
mm,  tg  =  0.8  mm,  L\  —  0.25  mm,  Le  =  100  mm,  Tv  =  300  K). 

A  comparison  with  the  numerical  data  reported  by  Stephan  and  Busse 
(1992)  has  also  been  made  for  ammonia  with:  Tv  —  300  K,  kw  =  221  W/(m2-K),  A' 
=  2  x  10~21  J,  a  =  1,  Li  =  0.5  x  10-5  m,  W  =  0.5  x  10~3  m,  7  =  45°,  tg  =  0.5  x  10-3 
m,  Rv  =  1  m,  R0  =  1.0015  m,  6j  =  $mtn  =  19.7°,  AT  =  1.31  K.  The  results  of  the 
comparison  are  listed  in  Table  3.2,  which  was  prepared  by  Stephan  and  Busse  (1992) 
except  for  the  data  of  the  present  analysis.  Tsv  is  the  temperature  of  the  vapor  side 
of  the  interface  and  Q1^ c  is  the  heat  flow  rate  per  unit  groove  length  in  the  region 
0  <  x  <  1  //m.  The  value  of  the  heat  transfer  coefficient  by  Stephan  and  Busse 
(1992)  was  23,000  W/(m2-K),  while  the  result  of  the  present  numerical  analysis  is 
17,385  W/(m2-K)  for  a  rough  surface  (for  f?r  =  0.02  ^um)  and  23,900  W/(m2-K)  for 
a  smooth  surface,  which  proves  the  validity  of  the  present  analysis. 
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Table  3.2:  Comparison  of  the  results  with  simplified  models 


Investigators 

Present 

Stephan 

and  Busse 

Schneider 

Shekriladze 

(1992) 

(1976) 

(1987) 

Assumption: 

Ts  >  Tv 

Ts  >  Tv 

TSv  >  Tv 

II 

S> 

II 

Ts  =  Tv 

Surface: 

rough 

smooth 

smooth 

smooth 

smooth 

smooth 

Tw-Tv  (K) 

1.31 

1.31 

1.31 

1.31 

- 

- 

he  (  W/cm2-K) 

1.74 

2.39 

2.3 

7.9 

6.9 

3.9 

qe  (W/cm2) 

2.56 

3.69 

3.0 

10.4 

- 

- 

QUQ't ot  (%) 

37 

38 

45 

94 

- 

- 
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The  influence  of  the  meniscus  contact  angle  on  the  local  heat  transfer  coeffi¬ 
cient  in  the  condenser  (configuration  by  Schlitt  et  al.  (1974))  is  demonstrated  in  Fig. 
3.7.  The  results  agree  qualitatively  with  those  obtained  by  Babenko  et  al.  (1981). 
The  increase  of  the  liquid  surface  curvature  causes  the  strong  decrease  of  the  heat 
transfer  coefficient,  where  a  sharp  maximum  occurs  in  the  vicinity  of  maximum  9men. 
In  this  location,  the  heat  transfer  coefficient  is  also  dependent  on  the  temperature 
drop  AT.  In  the  numerical  experiments  the  liquid  film  thickness  was  comparatively 
large  (Fig.  3.7(a))  and  the  interfacial  thermal  resistance  was  negligible  in  comparison 
with  that  of  the  film.  The  values  of  the  heat  transfer  coefficient  in  the  condenser 
based  on  the  outer  tube  diameter  for  the  ammonia  (Tv  =  250  K)  and  ethane  (Tv  = 
200  K)  heat  pipes  reported  by  Schlitt  et  al.  (1974)  are  7600  and  3300  W/(m2-K), 
respectively.  The  numerical  predictions  were  of  the  same  order  of  magnitude  as  that 
reported  by  Schlitt  et  al.  (1974). 

3.9  Conclusions 

The  results  of  the  numerical  simulation  of  heat  transfer  during  evaporation  and  con¬ 
densation  on  the  grooved  surfaces  of  heat  pipes  are  summarized  as  follows: 

1.  The  validity  of  the  present  mathematical  model  for  grooved  evaporators 
has  been  confirmed,  in  general,  by  the  comparisons  with  experimental 
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Figure  3.7:  Effect  of  the  meniscus  contact  angle  in  the  heat  pipe  condenser  on  (a) 
liquid  film  thickness  variation  along  the  surface  of  the  fin  top  (ammonia,  AT  =  1  K, 
Tv  =  250  K),  (b)  local  heat  transfer  coefficient  for  ammonia  (Tv  =  250  K),  (c)  local 
heat  transfer  coefficient  for  ethane  ( Tv  —  200  K) 
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data  by  Schlitt  et  al.  (1974),  Ivanovskii  et  al.  (1984)  and  the  numerical 
results  by  Stephan  and  Busse  (1992).  However,  more  detailed  information 
concerning  the  values  of  the  accommodation  coefficients  and  dispersion 
constants  is  needed. 

2.  Accounting  for  the  roughness  of  the  solid  surface  in  the  thin  evaporating 
film  region  resulted  in  a  decrease  of  the  heat  transfer  coefficient  by  up  to 
about  30%  in  comparison  to  that  obtained  for  a  smooth  surface  for  the  case 
when  the  accommodation  coefficient  was  set  equal  to  unity.  For  a  <C  1 
the  influence  of  the  surface  roughness  on  the  evaporative  heat  transfer 
coefficient  was  insignificant. 

3.  The  simplified  model  of  evaporative  heat  transfer,  were  it  was  assumed  that 
the  free  film  surface  curvature  in  the  microfilm  region  was  equal  to  that  in 
the  meniscus  region,  predicted  values  of  the  heat  transfer  coefficient  only 
up  to  5%  smaller  in  comparison  to  the  case  where  the  curvature  variation 
along  the  film  was  taken  into  account  (for  a  =  1). 

4.  The  value  of  the  local  evaporative  heat  transfer  coefficient  (for  a  fixed  9men) 
was  practically  independent  of  the  heat  flux  on  the  evaporator  external 
wall. 

5.  The  interfacial  resistance  significantly  influenced  the  value  of  the  evapo¬ 
rative  heat  transfer  coefficient.  Therefore,  the  more  advanced  expressions 
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for  this  resistance  during  high  intensive  evaporation  are  needed. 


6.  The  interfacial  resistance  and  the  disjoining  pressure  had  negligible  effects 
on  the  heat  transfer  in  grooved  condensers  because  of  the  large  thickness 
of  the  liquid  film,  while  the  effect  of  the  meniscus  contact  angle  was  sig¬ 
nificant. 

7.  Due  to  the  fact  that  the  greatest  part  of  the  heat  flow  through  the  liquid 
is  transferred  in  the  meniscus  region  where  the  liquid  film  is  comparatively 
small,  grooved  evaporators  are  capable  of  withstanding  comparatively  high 
heat  fluxes.  While  the  model  gives  the  values  of  the  temperature  drop  in 
the  liquid  during  evaporation,  it  provides  an  opportunity  to  predict  the 
onset  of  the  nucleate  boiling  in  grooved  evaporators. 

In  order  to  make  the  comparisons  with  experimental  data  more  profound, 
the  longitudinal  variation  of  the  meniscus  angle  should  be  taken  into  account.  That 
means  that  the  fluid  circulation  in  an  AGHP  should  be  considered  in  conjunction 
with  the  present  analysis. 
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Chapter  4 


ANALYSIS  OF  MINIATURE 
GROOVED  HEAT  PIPES 

4.1  Summary 

A  detailed  mathematical  model  of  axially-grooved  heat  pipes  (AGHP)  is  developed 
in  which  the  fluid  circulation  is  considered  along  with  the  heat  and  mass  transfer 
processes  during  evaporation  and  condensation.  The  predicted  results  obtained  are 
compared  to  existing  experimental  data.  Both  capillary  and  boiling  limitations  are 
found  to  be  important  for  the  flat  miniature  copper-water  heat  pipe,  which  is  capable 
of  withstanding  heat  fluxes  on  the  order  of  40  W/cm2  applied  to  the  evaporator  wall 
in  the  vertical  position.  The  influence  of  the  geometry  of  the  grooved  surface  on  the 
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maximum  heat  transfer  capacity  of  the  miniature  AGHP  is  demonstrated. 

4.2  Introduction 

Some  experimental  data  on  the  performance  characteristics  of  the  flat  water-copper 
heat  pipe  with  external  dimensions  2  x  7  x  120  mm  has  been  reported  by  Plesch  et 
al.  (1991).  The  present  chapter  deals  with  the  numerical  prediction  of  the  thermal 
characteristics  of  miniature  and  conventional  axially-grooved  heat  pipes  including 
the  maximum  heat  flux  on  the  outer  evaporator  wall.  This  heat  flux  is  restricted 
by  the  fluid  transport  limit  and  boiling  limitation.  The  superheat  of  the  evaporat¬ 
ing  liquid,  which  is  critical  for  the  onset  of  nucleate  boiling,  depends  on  the  heat 
transfer  coefficient  in  the  evaporator.  Therefore,  a  detailed  mathematical  model  is 
developed,  including  both  the  heat  transfer  through  liquid  films  and  axial  fluid  trans¬ 
port  phenomena,  which  provides  an  opportunity  to  determine  the  maximum  heat 
transfer  capacity  and  thermal  resistance  of  low-temperature  AGHPs.  The  present 
model,  which  incorporates  several  one-dimensional  boundary- value  problems,  has  the 
following  comparatively  new  features: 

•  The  heat  transfer  through  the  liquid  films  and  the  fin  between  grooves  (Fig.  4.1) 
in  both  the  evaporator  and  condenser  is  accounted  for  in  the  model,  which  is 
described  with  respect  to  the  disjoining  pressure,  interfacial  thermal  resistance 
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Figure  4.1:  Flat  miniature  axially-grooved  heat  pipe  cross  sections 


and  surface  roughness  or  curvature  in  the  same  manner  as  in  Chapter  3. 


•  The  effect  of  shear  stresses  at  the  free  liquid  surface  in  a  groove  due  to  the 
frictional  liquid- vapor  interaction  on  the  liquid  flow  is  taken  into  consideration. 

•  The  boiling  limitation  and  the  thermal  resistances  of  the  evaporator  and  con¬ 
denser  are  determined  by  accounting  for  the  longitudinal  distribution  of  the 
meniscus  curvature,  which  is  dependent  on  heat  load  and  heat  pipe  inclination 
angle. 

4.3  AGHP  Mathematical  Model 

The  mathematical  model  of  a  low-temperature  AGHP  is  developed  in  order  to  op¬ 
timize  the  performance  of  heat  pipes  and  grooved  evaporators  and  condensers  for 
cooling  systems.  The  proposed  model  is  a  significant  contribution  over  the  previ¬ 
ous  investigators’  attempts  (Kamotani,  1976b;  and  Vasiliev  et  al.,  1981),  where  the 
heat  transfer  during  evaporation  in  the  microfilm  and  boiling  limitation  were  not 
considered. 

The  present  model  is  developed  for  rectangular  and  triangular  (trapezoidal) 
grooves  in  circular  or  flat  (rectangular)  tubes,  but  other  heat  pipes  and  grooves  con¬ 
figurations  can  also  be  described  using  the  same  general  approach.  Heat  transfer 
processes  in  the  heat  pipe  container  and  working  fluid  were  considered  to  be  one- 
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dimensional  in  the  radial  direction,  such  that  axial  heat  conduction  was  neglected.  It 
was  also  assumed  that  no  puddle  flow  occurred  in  the  AGHP,  no  part  of  the  condenser 
was  blocked  by  the  liquid  working  fluid,  the  liquid  was  distributed  uniformly  between 
the  grooves,  and  the  fluid  flow  along  the  axis  was  related  to  the  capillary  potential 
gradient,  as  described  by  the  main  radius  of  curvature  of  the  liquid  in  a  groove. 

The  numerical  results  were  obtained  using  an  iterative  mathematical  proce¬ 
dure  which  involved  the  following  problems: 

1.  Heat  transfer  in  the  evaporating  film  on  a  rough  surface. 

2.  Heat  transfer  in  the  condensate  film  on  the  fin  top  surface. 

3.  Heat  conduction  in  a  metallic  fin  and  liquid  meniscus. 

4.  Fluid  circulation  in  the  AGHP. 

5.  Prediction  of  the  onset  of  the  nucleate  boiling. 

Since  the  first  three  problems  are  presented  in  detail  in  Chapter  3,  their  description 
is  omitted  here.  However,  these  interconnected  problems  are  included  in  the  math¬ 
ematical  procedure  along  with  the  fluid  circulation  and  boiling  limitation  equations 
presented  in  the  following  sections. 
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4.4  Fluid  Circulation  in  an  AGHP  and  the  Cap¬ 
illary  Limitation 

The  conservation  equations  for  the  steady-state  operation  of  an  AGHP  are  the  con¬ 
tinuity,  momentum  and  energy  equations  for  the  liquid,  vapor  and  wall,  and  the 
Laplace- Young  equation  for  the  radius  of  curvature  at  the  liquid- vapor  interface.  At 
any  axial  location,  the  following  mass  conservation  equation  must  hold  over  the  cross 
section  of  the  AGHP: 


WypyAy  =  Nw£P(A(  =  Nrnt  (4.1) 

where  rh£  is  the  mass  flow  rate  through  a  groove  and  N  is  the  number  of  grooves. 
The  average  liquid  and  vapor  velocities  in  the  axial  direction  in  equation  (4.1)  are 


Wi  =  J  Ja  w£(x,y)  dx  dy,  Wv  =  J  JA  wv(x>y)  dx  dy 


(4.2) 


An  AGHP  contains  a  definite  amount  of  a  working  fluid  Mt,  which  is  dis¬ 
tributed  in  accordance  to  the  following  relation,  provided  there  is  no  excess  liquid: 


Mt  —  N  j  p£A£  dz  +  J  pyAv  dz  -f  Ms  (4-3) 

where  Ms  is  the  mass  of  the  liquid  film  in  the  condenser  section  on  top  of  the  fins. 


rLt 
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The  conservation  of  axial  momentum  equation  for  an  incompressible  vapor 
flow  using  the  one-dimensional  boundary-layer  approximation  is 

d_ 

dz 

where  0V  is  the  momentum  flux  coefficient,  fv  is  the  friction  coefficient,  and  Dh,v 
is  the  hydraulic  diameter  of  the  vapor  channel.  These  coefficients  for  a  circular  va¬ 
por  channel  can  be  determined  using  the  results  of  previous  investigators  for  two- 
dimensional.  laminar  boundary-layer  flow  with  suction  or  injection.  For  the  range  of 
radial  Reynolds  numbers  usually  seen  in  miniature  AGHPs  (0  <  Rer  <  20)  in  the 
evaporator  and  adiabatic  sections,  the  results  given  by  Bankston  and  Smith  (1973) 
can  be  approximated  as  follows: 


(. Pv  +  Pv9z  sin  (p  -j-  pvfiywl )  =  -fv 


ZpvW2v 

■Dh^v 


(4.4) 


(/  Re)„  =  16  +  0.46Rer  —  0.017Re^,  f3v  =  1.33  —  0.005R.er 
In  the  condenser  section  (Bowman  and  Hitchcock,  1988;  Jang  et  ah,  1991) 


(/  Re)„  =  16(1.2337  -  0.2337exp(-0.0363Rer)][exp(1.2Ma)],  f3v  =  1.33 

The  velocity  of  the  liquid  phase  is  very  small  in  comparison  to  that  of  the  vapor 
flow,  so  the  interfacial  shear  for  the  vapor  is  computed  by  assuming  the  liquid  to  be 
stationary  (Longtin  et  ah,  1992).  For  a  flat  (rectangular)  vapor  channel  configuration 
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the  values  of  (/  Re)^  can  be  defined  using  the  following  equation  (Shah  and  Bhatti, 
1987) 


(/  Re)„  =  24(1  -  1.3553C  +  1.9467C2  -  1.7012C3  +  0.9564C4  -  0.2537C5) 

where  C  =  [t  —  2 [tw  +  tg)]/[N(W  4-  Xi)]  and  the  Reynolds  number  is  based  on  the 
hydraulic  diameter. 

The  axial  transport  of  condensate  in  an  AGHP  takes  place  in  the  grooves, 
where  most  of  the  liquid  resides.  The  thickness  of  the  liquid  film  on  the  fin  tops 
in  the  condenser  section  is  assumed  to  be  too  small  to  contribute  to  the  axial  mass 
transport.  Since  the  axial  Reynolds  number  for  liquid  flow  in  an  AGHP  is  usually 
less  than  10,  it  can  be  considered  to  be  viscous  and  quasi-one-dimensional.  Thus, 
the  conservation  of  momentum  equation  for  the  liquid  flow  in  a  groove  with  cross 
sectional  area  At  is 


7772 


dpe  .  2  ptwt 

—  +  ptg  sin  =  ft— - 

dz  Dh,£ 


(4.5) 


The  values  of  ft  can  be  taken  from  the  numerical  solution  of  the  two-dimensional 
Stokes  equation  with  boundary  conditions  which  take  into  account  shear  stress  in 


the  liquid  at  the  interface  due  to  liquid- vapor  frictional  interaction  (Vasiliev  et  ah, 
1981;  Schneider  and  DeVos,  1980).  The  expression  for  the  friction  factor,  given  by 
Schneider  and  DeVos  (1980)  for  rectangular  grooves  is  used  in  the  presented  model 
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because  of  its  reliability  and  accuracy.  For  the  case  of  a  planar  free  liquid  surface 
geometry 


d_ 

dz 


(pewtAe) 


-h0,e{T0,e-Tv)(l  +  ^\  , 


h0,e  A 


0,  (adiabatic  section) 

Nhfgko’ci'To’c~  Tv){l  +  X)  ’ 


0  <z  <  Le 
Le  z  Le  d-  L a 

Le  +  La  <  z  <  Lt 


(4.6) 


he(z)  and  hc(z)  are  the  mean  local  effective  heat  transfer  coefficients  between  the 
external  surface  of  the  casing  and  vapor  through  the  wall  and  liquid  films  (Khrustalev 
and  Faghri,  1993).  h0<e  and  h0<c  are  the  heat  transfer  coefficients  between  the  wall 
and  the  ambient.  The  vapor  temperature  is  denoted  by  Tv.  For  a  given  axial  heat 
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load  function  Q(z),  equation  (4.6)  can  be  rewritten  as 


d  _  1  dQ(z) 

tlz  PtU’1  1  -  hj,N  dz 


(4.7) 


The  interfacial  radius  of  the  meniscus  curvature  is  related  to  the  pressure 
difference  between  the  liquid  and  vapor  by  the  Laplace- Young  equation,  which,  in 
differential  form,  is 


dpt  dpv  d  /  a 
dz  dz  dz  \Rm 


(4.8) 


The  mean  axial  liquid  and  vapor  velocities  at  the  evaporator  and  condenser 
end  caps  are 


Wi\z=o  =  wv\z=0  =  0  (4.9) 

m\ Z-Lt  =  wv\z=Lt  =  0  (4.10) 

which  are  imposed  by  the  axial  heat  load  function  Q(z).  The  vapor  and  liquid  pres¬ 
sures  at  the  evaporator  end  cap  are 


Pv  \z~0  —  Pv 0? 


PtlzzzO  Pv  0 


(4.11) 


The  values  of  pv o  and  Rm o  are  to  be  determined  using  additional  conditions.  When 
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solving  equation  (4.6),  where  Tv(z)  is  also  unknown  and  depends  on  pV:  the  value  of 
pv o  must  be  chosen  such  as  to  satisfy  equation  (4.10).  When  Q(z)  and  Tv0  are  known, 
Pvo  =  Pv{Tv o),  which  can  be  obtained  using  the  saturation  tables.  When  there  is  an 
insufficient  liquid  fill  volume,  the  radius  of  curvature  of  the  meniscus  at  the  evaporator 
end  cap,  Rm 0,  is  found  using  equation  (4.3),  since  the  value  of  Rm0  influences  the  axial 
distribution  of  liquid.  In  general,  the  radius  of  curvature  of  the  meniscus  is  bounded 
by  Rm, min  <  Rm o  <  Rv  Since  in  the  presented  numerical  experiments  the  situation  of 
a  sufficient  liquid  charge  was  considered,  the  following  conditions  for  Rm0  definition 
were  used 


min{  }  e>  >  rnin{  D  ^  }  0 


(4.12) 


Rm(z)J  Rv  ’  “~~lRm(z)- 

for  circular  and  flat  heat  pipes,  respectively.  Equation  (4.12)  implies  that  for  a  given 
Q(z),  the  value  of  Rm o  which  provides  the  equality  of  Rm  to  the  vapor  channel  radius 
(for  the  circular  tube  case)  at  a  single  location  on  z  along  the  AGHP  should  be  chosen, 
while  equations  (4.4)-(4.11)  are  satisfied  along  the  entire  length  of  the  AGHP.  Note 
that  this  point  is  not  always  situated  at  the  condenser  end  cap. 

The  solution  of  equations  (4.4)-(4.11)  gives  the  longitudinal  distribution  of 
the  meniscus  radius  Rm(z)  and  hence  the  meniscus  contact  angle  #men(2),  which  are 
related  for  contact  angles  of  #men  >  ^min  by 
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(4.13) 


^men(^)  =  arCCOS 


W 


L  Rm(z) 


-7 


For  rectangular  grooves  Rm>m m  is  defined  by  equation  (4.13)  provided  0men  =  O^n, 
where  6^ n  is  the  minimum  wetting  contact  angle,  which  is  fixed  for  a  specific  working- 
fluid/container  combination  (Stepanov  et  ah,  1977).  This  value  of  Rm> mjn  along  with 
equation  (4.12)  determines  the  maximum  heat  transfer  through  an  AGHP  (capillary 
limit)  at  which  equations  (4.4)-(4.11)  and  the  boundary  conditions  are  still  satisfied. 
For  triangular  (trapezoidal)  grooves  in  the  situation  9men  =  ^min.  equation  (4.13)  gives 
the  relation  between  the  meniscus  radius  and  the  width  of  the  liquid  cross  section. 

When  the  longitudinal  vapor  pressure  distribution,  pv(z),  and  the  tempera¬ 
ture  at  the  evaporator  end  cap,  Tv o  (for  the  case  of  an  unknown  heat  load  function) 
are  found,  the  corresponding  temperature  distribution  along  the  vapor  flow,  Tv(z), 
can  be  calculated  using  the  Clausius-Clapeyron  equation. 


4.5  Boiling  Limitation  and  Heat  Transfer  Coeffi¬ 
cients 

Nucleation  within  the  wick  is  undesirable  for  wicked  heat  pipe  operation  because  the 
bubbles  can  obstruct  the  liquid  circulation  and  hence  cause  hot  spots  on  the  evapo¬ 
rator  walls.  An  estimation  of  the  vapor  bubble  departure  diameter  during  nucleate 
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boiling  of  water  was  made  using  the  correlations  obtained  by  Jensen  and  Memmel 
(1986).  The  bubble  departure  diameter  was  much  larger  than  the  groove  width  of 
the  considered  water-copper  heat  pipes,  which  means  that  the  normal  operation  of  an 
AGHP  in  the  situation  when  nucleate  boiling  occurs  in  the  evaporator  is  not  realistic 
except  for  some  special  cases.  Thus,  for  heat  pipe  performance  optimization,  it  is 
necessary  to  predict  the  value  of  the  heat  flux  which  initiates  vapor  bubble  formation 
in  the  working  liquid.  This  task  can  be  subdivided  into  two  independent  parts: 

1.  Estimation  of  the  superheat  of  the  liquid  in  the  wick  (Tb0t,e  —  Tv)  which 
causes  generation  of  vapor  bubbles  at  the  wall-wick  interface.  Here  Tb0t,i 
corresponds  to  the  plane  x  =  tg  (Fig.  4.1(b)). 

2.  Computation  of  the  heat  transfer  coefficient  in  the  evaporator  (he)  for  the 
situation  when  no  boiling  occurs,  which  relates  the  values  of  the  heat  flux 
(■ qe )  and  temperature  drop  across  the  wick  (That/  —  Tv). 

Critical  liquid  superheat.  A  review  on  heterogeneous  nucleation  and 
bubble  growth  in  liquids  has  been  recently  given  by  Carey  (1992).  Since  the  inner 
surface  of  heat  pipes  are  invariably  covered  with  machine- formed  pits,  cavities  and 
scratches,  vaporization  usually  begins  on  these  irregularities.  The  process  of  bubble 
formation  depends  upon  the  solid  surface  microrelief,  wetting  contact  angles  and 
liquid  properties.  While  theoretical  prediction  of  the  superheat  which  causes  boiling 
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is  extremely  difficult,  some  empirical  correlations  should  be  involved  in  the  analysis. 
Note  that  the  vapor  bubble  embryo  diameter  is  much  smaller  than  the  groove  width, 
therefore,  some  of  the  correlations  obtained  for  planar  surfaces  can  be  used  in  the 
present  analysis. 

Lorenz  et  al.  (1974)  reported  experimental  data  for  boiling  of  four  different 
liquids  on  a  copper  surface  with  a  #240  (sandpaper)  finish.  In  these  experiments,  the 
number  of  active  sites  per  square  centimeter  (n),  which  increased  rapidly  with  super¬ 
heat,  were  determined  by  visual  counting  at  different  superheat  values.  Experimental 
results  by  Lorenz  et  al.  (1974)  can  be  approximated  by  the  following  equation 


h  =  -  C2 


bot,^  Ty)py  111 


(4.14) 


2crTv  '  2 Rn 

where  for  water  C\  =  77,  C2  —  19xl06  and  1  <  n  <  6,  and  for  R-113,  methanol 
and  benzene  C\  =  39.6,  C2  =  19.6x10s  and  0.5  <  n  <  10.  Note  that  in  the  analysis 
by  Lorenz  et  al.  (1974)  Ts at(p^)  was  used  instead  of  Tv  and  l/Rm  was  equal  to  zero. 
The  following  relation  between  Tsa.t(pi)  and  Tv  was  used  in  the  derivation  of  equation 
(4.14)  (Carey,  1992,  p.181) 


T, ,  -  TUn)  =  t  ~~~r> 

^  /  g  Pv-^m 

with  the  assumptions  that  pv  <C  pt  and  0 / piRm  1-  For  a  known  superheat 

(Tbot/  —  Tv),  equation  (4.14)  indicates  whether  boiling  occurs  at  a  given  point  along 
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the  evaporator  length,  which  is  characterized  by  the  local  radius  of  curvature  of  the 
meniscus,  Rm.  The  values  of  the  heat  flux  on  the  external  surface  of  the  wall  and  the 
superheat  are  related  as  follows 

qe  =  he,bot(TboM  -  Tv)  (Rv+£  +  t')  (4-15) 

where  /ie,bot  is  the  local  effective  heat  transfer  coefficient  between  the  saturated  vapor 
and  the  bottom  of  a  groove,  which  is  generally  dependent  on  the  value  of  Rm  at  this 
location. 

Heat  transfer  coefficients.  The  coefficients  of  heat  transfer  during  evap¬ 
oration  and  condensation  on  the  grooved  surface  were  calculated  using  the  mathe¬ 
matical  model  developed  in  Chapter  3.  It  is  essential  to  mention  that  in  the  model 
given  in  Chapter  3  the  thin  film  formation  on  the  fin  surface  for  both  evaporation 
and  condensation  was  described  by  the  following  equation: 


ke(Tw  -  Ts) 


(4.16) 


hfg  Pt$ 

where  s  is  the  coordinate  along  the  solid-liquid  interface  where  the  temperature  is  Tw, 
and  6  is  the  local  film  thickness.  The  temperature  of  the  liquid-vapor  interface  T$  is 
affected  by  the  disjoining  and  capillary  pressures,  and  also  depends  on  the  value  of 
the  interfacial  resistance,  which  is  defined  for  the  case  of  a  comparatively  small  heat 
flux  by  the  following  relation 
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2a  \  ^Jg  [  Pv  (.Psat) 


2-<*/  J2^Rgi^/Tv  VTs 


(4.17) 


pv  and  (psat)^  are  the  saturation  pressures  corresponding  to  Tv  in  the  bulk  vapor  and 
at  the  thin  liquid  film  interface,  respectively.  The  extended  Kelvin  equation  which 
relates  p<*,  psa.t(Ts )  (normal  saturation  pressure  corresponding  to  T$),  and  (psa t)s  was 
also  used  in  the  model.  The  following  equations  for  the  disjoining  pressure  were  used 
in  the  present  model  for  nonpolar  liquids  and  for  water,  respectively: 


Pd  =  -A' S' 


(4.18) 


Pd  =  f>iR,Ts  In  o  I  j-, 


(4.19) 


where  a  =  1.534  and  b  =  0.0243.  Solution  of  the  heat  conduction  equation  in  the  fin 
and  liquid  in  conjunction  with  the  equations  for  the  thin  film  formation  for  a  given 
qe  gives  the  temperature  of  the  fin  bottom  as  T\x=tg  —  Tb0t/-  The  local  heat  transfer 
coefficient  from  the  bottom  of  the  groove  to  the  vapor  is  defined  as 


he,bot  — 


Qe  (  Ro 

[T|x=ts  —  Tv]  yiik  +  tt  +  tg 


(4.20) 


The  local  heat  transfer  coefficient  from  the  external  surface  of  the  heat  pipe  to  the 


vapor  is 
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-1 


(4.21) 


h  =  A  ln  f  ^  )  |  _ 1  (  R”  ) 

lkw  \Ry  +  tt  +  tg  )  he  ,bot  \RV  +  tt  +  tg  J 
where  the  thermal  resistance  of  a  circular  tube  wall  is  taken  into  account.  It  was  found 
in  the  results  by  Khrustalev  and  Faghri  (1993)  that  for  a  =  1  the  values  of  he  for  the 
rough  solid  surface  were  up  to  30%  less  than  for  the  smooth  solid  surface.  However, 
for  the  case  of  small  a  (for  example,  a.  =  0.05)  the  influence  of  the  roughness  on  the 
evaporative  heat  transfer  coefficient  was  insignificant.  The  characteristic  roughness 
size  Rr  in  the  referred  results  varied  from  0.02  to  1  ft m.  Obviously,  equations  (4.15), 
(4.20)  and  (4.21)  are  valid  also  for  a  planar  evaporator  in  which  R0  — *•  oo,  Rv  — >  co 
and  R0  =  Rv  +  +  tg  +  tw. 

In  the  present  analysis  the  values  of  he,bot  were  calculated  for  every  point 
along  the  axial  direction  in  the  evaporator  of  the  modeled  AGHPs  for  a  given  total 
heat  load  Qa ,  and  equation  (4.14)  was  used  with  the  corresponding  Rm  to  predict  the 
onset  of  nucleate  boiling  (n  >  1).  The  mean  temperature  drop  in  the  evaporator  was 
calculated  using  the  following  expression: 
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4.6  Numerical  Treatment 


For  the  fluid  circulation  in  the  AGHP  (equations  (4.4)-(4.13)),  the  four  main  variables 
dependent  on  2  to  be  found  are 

K  =  i?"1:  curvature  of  the  free  surface  of  the  meniscus, 
rh  =  wtpeAf.  mass  flow  rate  of  liquid  in  one  groove, 
pp.  liquid  pressure, 
pv:  vapor  pressure. 

The  system  of  four  ordinary  differential  equations,  equations  (4.4),  (4.5),  (4.7)  and 
(4.8),  was  solved  with  the  respective  boundary  conditions  using  the  fourth-order 
Runge-Kutta  procedure  and  the  shooting  method  (on  parameter  Rm 0).  For  each 
of  the  four  functions  the  relative  error  was  less  than  0.0001%.  While  in  the  present 
paper  the  AGHP  for  the  case  of  a  given  axial  heat  load  function  Q{z)  was  modeled, 
the  structure  of  the  numerical  AGHP  simulation  was  as  follows: 

1.  The  fluid  circulation  along  the  AGHP  was  solved  with  the  known  function 
Q(z)  and  Tv 0,  which  gave  Rm(z),  0men (z)  and  Tv(z). 

2.  For  every  point  on  z,  the  heat  transfer  in  the  evaporator  and  condenser 
thin  film  regions  along  with  the  heat  conduction  in  the  fin  were  solved  and 
the  functions  he{z),  he^ol{z),  hc(z)  and  TWi0(z )  were  determined. 

3.  The  number  of  active  sites  n  was  calculated  for  every  point  on  z  in  the 
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evaporator  with  the  corresponding  Rm  to  determine  whether  boiling  oc¬ 
curred. 

4.  The  values  of  the  temperature  drops  in  the  evaporator  and  condenser  were 
found. 

5.  Steps  1-4  were  repeated  several  times  for  larger  values  of  Q a  until  the 
capillary  or  boiling  limitation  was  indicated  (with  the  relative  error  less 
than  1%). 

4.7  Results  and  Discussion 

To  verify  the  numerical  results  obtained,  the  experimental  data  provided  by  Schlitt  et 
al.  (1974)  were  first  used.  Therefore,  part  of  the  presented  results  refer  to  the  AGHP 
with  the  following  geometry:  Lt  =  0.914  m,  Lc  =  0.152  m,  0.15  <  Le  <  0.343  m,  W 
=  0.305  mm,  tg  =  1.02  mm,  -  0.215  mm,  Rv  =  4.43  mm,  R0  -  7.95  mm,  7  =  0, 
N  =  27,  tt  =  0.  The  working  fluids  were  ammonia  and  ethane,  the  casing  material 
thermal  conductivity  was  assumed  to  be  kw  =  170  W/(m-K),  and  dispersion  constant 
A'  =  10-21  J.  The  axial  heat  distribution  was  specified  as 
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Q={ 


(4.23) 


Q  a.Z  /  Le-, 

0  <  z  <  Le 

Qa, 

Le  Z  <  Le  ~t~  L^i 

*)’ 

Le  T  La  5;  z  <  Lt 

and  equation  (4.7),  instead  of  equation  (4.6),  was  used  in  the  numerical  procedure. 
The  longitudinal  distribution  of  the  meniscus  contact  angle,  which  is  influenced  by 
the  inclination  angle  and  the  fluid  pressure  variation  (Ap),  is  shown  in  Fig.  4.2(a) 
for  ammonia  (Le  =  0.343  m,  Tv  =  250  K).  For  positive  values  of  the  inclination  angle 
V5  (when  the  condenser  end  is  elevated)  the  points  of  minimum  and  maximum  liquid 
surface  curvature  (the  so-called  “dry”  and  “wet”  points)  were  shifted  from  the  ends  of 
the  heat  pipe  towards  the  adiabatic  section.  As  a  result,  the  meniscus  contact  angle 
distributions  in  the  heat  loaded  sections  were  even  more  uniform  than  those  in  the 
horizontal  case.  Note  that  the  meniscus  angles  in  the  condenser  differ  very  slightly 
from  the  maximum  value,  and  their  values  in  the  evaporator  can  be  almost  the  same 
as  that  in  the  condenser  for  positive  p.  The  corresponding  longitudinal  distributions 
of  the  local  heat  transfer  coefficients  in  the  evaporator  and  condenser  and  also  wall 
and  vapor  temperature  variations  are  shown  in  Figs.  4.2(b)  and  4.2(c).  The  local  heat 
transfer  coefficient  in  the  middle  of  the  condenser  was  about  two  times  larger  than 
at  the  entrance  or  end  cap,  which  resulted  in  the  external  wall  temperature  variation 
shown  in  Fig.  4.2(c).  In  the  evaporator  section,  the  variation  of  the  local  heat  transfer 
coefficient  was  weaker  (for  heat  loads  which  are  not  close  to  the  maximum),  so  the 
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Figure  4.2:  Performance  characteristics  of  the  ammonia-Al  heat  pipe  (Tv  =  250  K) 
(a)  meniscus  contact  angle  and  fluid  pressure,  (b)  local  heat  transfer  coefficients,  (c) 
wall  and  vapor  temperatures 
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wall  temperature  profile  was  very  smooth.  The  temperature  drop  along  the  vapor 
flow  was  less  than  0.01  K. 

The  comparison  of  the  computational  results  with  the  experimental  data 
on  heat  transfer  for  the  AGHPs  evaporators  by  Schlitt  et  al.  (1974)  and  for  the 
planar  water-copper  evaporator  by  Ivanovskii  et  al.  (1984)  is  presented  in  Fig.  4.3 
for  ammonia  ( Le  =  0.343  m,  Tv  =  250  K),  ethane  ( Le  =  0.15  m,  Tv  —  200  K)  and 
water  (W  =  0.17  mm,  tg  =  0.8  mm,  Li  =  0.25  mm,  Le  =  100  mm,  Tv  =  300  K; 
the  heated  surface  was  of  about  40  cm2).  The  agreement  of  the  results  is  good  for 
a  <C  1  since  it  was  noted  by  Carey  (1992,  p.  119)  that,  for  some  substances  (ethanol, 
methanol,  water,  etc.),  the  accommodation  coefficient  had  been  found  to  have  very 
small  values  (0.02  to  0.04)  in  the  experiments  of  Paul  (1962).  The  physical  reason 
for  low  values  of  a  in  the  microfilm  region  of  the  evaporator  can  be  the  concentra¬ 
tion  of  the  contaminants  which  usually  exist  in  a  heat  pipe  in  this  region.  While  the 
experimental  data  for  the  temperature  drops  for  ethane  were  still  higher  than  those 
predicted,  more  precise  values  of  the  accommodation  coefficient  and  dispersion  con¬ 
stants  are  needed.  For  the  copper-water  evaporator  the  predicted  boiling  limit  heat 
flux  was  37  W/cm2,  and  no  boiling  was  reported  by  Ivanovskii  et  al.  (1984),  even 
though  the  liquid  superheat  exceeded  20  K  at  qe  —  20  W/cm2  in  the  experiments. 
For  a  larger  number  of  grooves  per  unit  width  such  as  in  the  AGHP  by  Plesch  et 
al.  (1991),  the  heat  transfer  coefficient  should  be  higher  because  of  the  increase  in 
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Figure  4.3:  Temperature  drop  in  the  evaporator  versus  heat  load  (a)  ammonia-Al 
heat  pipe  (Tv  =  250  K),  (b)  ethane-Al  heat  pipe  ( Tv  —  200  K),  (c)  water-copper 
planar  evaporator  (Tv  =  300  K)  and  flat  miniature  AGHP  ( Tv  =  378  K) 
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the  thin  film  region  over  the  same  area.  The  characteristic  values  of  the  evaporative 
heat  transfer  coefficients  for  the  evaporator  by  Ivanovskii  et  al.  (1984)  were  10,000 
W/(m2-K),  while  for  the  miniature  AGHP  by  Plesch  et  al.  (1991)  the  predicted  val¬ 
ues  of  he  were  larger  than  40,000  W/(m2-K).  Thus,  the  boiling  limitation  should  take 
place  at  larger  heat  fluxes  for  narrower  grooves. 

The  agreement  between  the  predicted  temperature  drop  in  the  condenser 
and  the  experimental  results  by  Schlitt  et  al.  (1974)  can  be  seen  from  Fig.  4.4 
for  the  ammonia  (Tv  =  250  K)  and  ethane  (Tv  =  200  K)  heat  pipes.  As  the  heat 
load  increased,  the  thermal  resistance  of  the  condenser  increased  and  the  evaporator 
thermal  resistance  decreased  (see  Figs.  4.3  and  4.4)  due  to  the  changes  in  the  meniscus 
contact  angle  longitudinal  distribution. 

The  numerically  obtained  values  for  the  maximum  heat  transfer  capacity 
of  the  AGHPs  have  been  also  compared  with  the  experimental  data  by  Schlitt  et  al. 
(1974)  for  ammonia  (Fig.  4.5(a),  Le  =  0.3  m,  Tv  =  203  K)  and  ethane  (Fig.  4.5(b), 
Le  =  0.15  m,  Tv  =  200  K)  within  the  assumption  that  =  0.  Some  disagreement  of 
the  results  for  comparatively  large  inclination  angles  (elevation)  when  the  evaporator 
end  is  elevated  can  be  explained  by  the  complicated  meniscus  shape  at  the  evaporator 
end,  which  is  not  described  by  the  present  model.  However,  for  smaller  tilt  angles  the 
agreement  is  satisfactory. 

The  second  heat  pipe  configuration  considered  here  is  the  AGHP  investi- 
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Figure  4.4:  Temperature  drop  in  the  condenser  versus  heat  load:  (a)  ammonia-Al 
heat  pipe  (Tv  =  250  K),  (b)  ethane- A1  heat  pipe  (Tv  =  200  K) 
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Figure  4.5:  Maximum  heat  transfer  versus  elevation  height  (a)  ammonia-Al  heat  pipe 
{Tv  -  203  K),  (b)  ethane- Al  heat  pipe  (Tv  =  200  K),  (c)  water-copper  flat  miniature 
heat  pipe 
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gated  experimentally  by  Plesch  et  al.  (1991).  The  cross  sections  are  shown  in  Fig. 
4.1  and  the  AGHPs  geometrical  characteristics  are:  W  =  0.06  mm,  Li  =  0.02  mm,  tg 
=  0.24  mm,  t  =  2  mm,  tw  =  0.22  mm,  Le  =  20  mm,  La  -  80  mm,  Lc  =  20  mm,  N  = 
80.  The  casing  material  was  copper  and  the  working  fluid  was  pure  water.  Therefore, 
the  minimum  contact  angle  was  set  to  ^mjn  =  33°  (Stepanov  et  al.,  1977).  The  heaters 
were  situated  at  both  sides  of  the  evaporator  with  an  overall  heated  surface  of  2.56 
cm2.  The  numerical  results  for  the  maximum  heat  flow,  Qmax,  transferred  by  the  flat 
miniature  AGHP  are  shown  in  Fig.  4.5(c),  where  it  can  be  seen  that  Qmax  increased 
with  the  operating  temperature.  Flowever,  for  Tv  >  100°C,  it  was  restricted  by  the 
boiling  limitation.  To  validate  the  numerical  results,  a  comparison  with  two  of  the 
most  representative  series  of  the  experimental  data  by  Plesch  et  al.  (1991)  was  made, 
as  shown  in  Fig.  4.5(c).  It  should  be  noted  that  the  experimental  results  by  Plesch 
et  al.  (1991)  were  given  (and  plotted  in  Fig.  4.5(c))  versus  the  evaporator  surface 
temperature,  and  it  was  not  mentioned  directly  if  the  data  for  the  measured  heat  flow 
rate  were  at  a  maximum.  The  numerical  predictions  of  the  capillary  limit  for  both 
the  horizontal  (<p  =  0)  and  vertical  (<p  =  90°)  orientations  follow  the  experimental 
data.  The  temperature  drop  along  the  vapor  flow  was  less  than  0.4  K  in  the  numeri¬ 
cal  experiments  while  the  temperature  drop  in  the  evaporator  reached  several  degrees 
(Fig.  4.3(c)).  The  influence  of  the  grooved  surface  geometry  on  the  maximum  heat 
transfer  through  the  flat  AGHP  with  the  external  dimensions  2  x  7  x  120  mm  is 
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shown  in  Figs.  4.6  and  4.7  for  the  vertical  and  horizontal  orientations,  respectively. 
While  the  values  of  W  and  tg  were  changed,  the  fin  half- width  and  the  wall  thickness 
were  kept  constant  at  Lx  —  0.02  mm  and  tw  =  0.22  mm.  The  number  of  grooves 
N  was  different  for  different  W,  while  the  external  dimensions  were  the  same.  Also, 
the  vapor  channel  height  was  dependent  on  tg  because  of  the  same  reason.  For  the 
vertical  orientation  at  Tv  =  55°  C,  the  capillary  limit  restricted  the  maximum  heat 
flux  in  the  evaporator,  which  had  a  maximum  value  of  about  15  W/cm2  for  W  = 
0.03  to  0.04  mm  (Fig.  4.6(a)).  For  Tv  =  105°C  both  the  capillary  and  boiling  limits 
occurred  as  shown  in  Fig.  4.6(b),  and  the  optimal  groove  half  width  was  about  0.05 
mm  with  an  absolute  maximum  heat  flux  of  40  W/cm2.  In  the  horizontal  orientation 
the  maximum  heat  fluxes  were  8  W/cm2  (for  W  =  0.03  mm,  Tv  =  55°C)  and  25 
W/cm2  (for  W  =  0.04  mm,  Tv  =  105°C)  as  shown  in  Fig.  4.6. 

It  is  anticipated  that  in  order  to  reach  higher  heat  fluxes,  more  advanced 
heat  pipe  configurations  should  be  used. 

4.8  Conclusions 

The  numerical  results  of  the  mathematical  model  describing  the  fluid  flow  and  heat 
transfer  in  the  conventional  and  miniature  axially-grooved  heat  pipes  are  summarized 
as  follows: 
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Figure  4.6:  Maximum  heat  transfer  of  the  flat  miniature  heat  pipe  versus  groove 
depth  (vertical  orientation)  (a)  Tv  =  55°,  (b)  Tv  =  105° 


1.  The  validity  of  the  present  mathematical  model  for  axially- grooved  heat 
pipes  has  been  confirmed,  in  general,  by  comparisons  with  existing  exper¬ 
imental  data.  However,  more  detailed  information  concerning  the  values 
of  the  accommodation  coefficients  and  dispersion  constants  is  needed. 

2.  An  increase  of  the  heat  load  decreased  the  evaporator  thermal  resistance 
and  increased  the  condenser  thermal  resistance  due  to  changes  in  the  lon¬ 
gitudinal  meniscus  angle  distribution  along  the  AGHP  in  the  horizontal 
position. 

3.  The  maximum  heat  flux  in  the  evaporator  of  the  flat  miniature  copper- 
water  AGHP  was  restricted  by  both  the  capillary  and  boiling  limitations 
in  the  interval  of  the  operating  temperature  from  Tv  —  90  to  120°C. 

4.  Due  to  the  fact  that  the  greatest  part  of  the  heat  flow  through  the  liquid  is 
transferred  in  the  region  where  the  liquid  film  thickness  is  comparatively 
small,  grooved  evaporators  are  capable  of  withstanding  high  heat  fluxes. 
The  optimized  copper-water  flat  AGHP  with  external  dimensions  2x7 
x  120  mm  could  operate  at  maximum  heat  fluxes  in  the  evaporator  of  25 
and  40  W/cm2  for  the  horizontal  and  vertical  orientations,  respectively, 
provided  that  the  operating  temperature  was  high  enough  (  Tv  =  90  to 
120°C). 
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